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1 Introduction 

We present a method for calculating the zeta function of a smooth projec- 
tive variety over a finite field which proceeds by induction on the dimension. 
Specifically, we outline an algorithm which reduces the problem of calculating a 
numerical approximation for the action of Frobenius on the middle-dimensional 
rigid cohomology of a smooth variety, to that of performing the same calculation 
for a smooth hyperplane section. We present in detail the main new algorith- 
mic ingredient under some simplifying assumptions, and give full details of our 
algorithm for calculating zeta functions for some specific surfaces; we call it the 
"fibration algorithm". We have implemented the fibration algorithm for these 
surfaces over prime fields using the Magma programming language, and present 
some explicit examples which we have computed. 

To illustrate the main idea behind our approach, we begin by outlining 
the proof given by Dcligne of the Riemann hypothesis for a smooth projective 
variety X over the finite field Specifically, the statement that for each 

< * < 2 dim(X) the action of the Frobenius endomorphism on the €-adic etale 
cohomology space Hl^{X, Q^) has eigenvalues of complex absolute value q^/'^. 

Let X C P be a smooth projective variety of dimension n -I- f > 1 defined 
over the finite field F^. Denote by P the dual projective space whose points t 
correspond to hyperplanes Ht in P, and let D he a, line in P. Let X G X x D 
denote the set of points {x,t) such that x € Ht- Projection on the first and 

second coordinates yields maps X ^ X D. The fibre of / at i e D is 
the hyperplane section Xt — X O Ht oi X. For sufficiently general D these 
maps define a Lefschetz pencil 9, (5.1)] (one may need to change the projective 
embedding first (5.7)]). 

The action of the Frobenius endomorphism on the ^-adic etale cohomology 
H*^(X, Q() may be studied via this Lefschetz pencil. In particular, assuming the 
result holds for smooth curves and arguing by induction on the dimension n+l, 
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one can reduce the proof of the Riemann hypothesis for X to the case of the 
Frobenius action on the middle-dimensional cohomology space H'^^^{X,Q^i). 
The Leray spectral sequence for / and further inductive arguments now reduce 
the proof of the Riemann hypothesis to the case of E^'^^ := H^^^D , R"- f^Q^) 
(7.1)]. Specifically, one must prove that the Frobenius acting on this finite 
dimensional Q^-vector space has eigenvalues which have complex absolute value 
q{n+i)/2^ This is the "core problem" and it requires considerable ingenuity. 

In this article, we are interested in computing the eigenvalues of Frobenius, 
rather than proving that they verify Weil's conjecture. However, it should be 
possible to bring to bear upon this computational problem the above geomet- 
ric machinery. Specifically, one expects that the geometric techniques which 
Deligne used in his reduction to the core problem can be made algorithmic. 
However, even once this is done, one is still faced with a difficult problem viz. 
calculation of the Frobenius action on -E2 "t- The present author has no idea 
on how this might be achieved. However, the sketch of Deligne's proof can be 
presented in the terms of rigid cohomology, rather than £-adic etale cohomol- 
ogy, and this theory is much more amenable to computation. In this article 
we present an algorithmic solution to the analogous "core problem", at least 
under certain simplifying assumptions. The principal novelty of this algorith- 
mic technique is that it proceeds by induction on the dimension. Specifically, 
the calculation of a matrix for the action of Frobenius on the rigid cohomolog- 
ical analogue i?2'"ig requires as input a matrix for the action of Frobenius on 
H"^g{Xt) for some hyperplane section Xt of X. So we can show that for the 
purposes of computation, the "core problem" of calculating Frobenius in the 
middle dimension can be efficiently reduced to that of a single instance of the 
problem one dimension lower down. In our method the base case of curves is 
handled using Kedlaya's algorithm [TSj . 

We note in passing that for smooth projective hypersurfaces (of odd di- 
mension) Deligne's solution of the "core problem" can be applied in a different 
manner, viz. rather than fibring the hypersurface X in a Lefschetz pencil, one 
can embed it as a fibre in such a pencil |H1 (5.12)],^^]. Such an approach to 
calculating zeta functions was taken by the author in the "deformation algo- 
rithm" |21| . From a computational point of view, this latter approach has the 
disadvantage that the "total space" under consideration has dimension one more 
than the hypersurface itself. This impacts somewhat on the complexity of the 
"deformation algorithm". Specifically, the time/space complexity in terms of 
the "middle betti number" dim H^!'^^ (X) is rather high. Our new approach, of 
fibring the original variety, though more complicated, does appear better from 
the point of view of complexity dependence on the "betti numbers" . 

The algorithm presented in this paper uses the main technique developed 
for the "deformation algorithm" , combined with a "higher rank" generalisation 
of Kedlaya's algorithm. Although our recursive approach was conceived as a 
general purpose algorithm, our implementation and complexity analysis for some 
surfaces suggest it is likely to be of most use for surfaces which can be fibered 
into low genus curves. Specifically, for the surfaces we consider in Sections 
13 IHl and ini if one fixes the genus g of the generic fibre of the fibration, then 
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the asymptotic complexity of our algorithm is quasi-quartic in the middle betti 
number, with quasi-cubic space requirement. In fact, the complexity in this 
case is comparable to that in the original algorithm of Kedlaya ^Sl j only in this 
case we have surfaces rather than curves (Theorem I8.6|l . The dependence on 
the genus g itself is roughly comparable to that in the "deformation algorithm" 
for curves; see the end of Section 

We now outline the contents of the various sections in this paper. In Section 
|21we define the zeta function of a variety and explain the computational problem 
which pertains to them. In Section |2| we give the main definitions from rigid 
cohomology which we shall need, and define the specific computational problem 
on which we shall focus (Problem viz., calculation of Frobenius on the 
space "io- Neither Section |2] nor contains any original contribution. 

Section considers an "abstract" version of the main computational prob- 
lem, and proves a number of theorems relevant to its solution (Theorems 14.21 
and I4.8|l . The main theorem stated in this section (Theorem 14. 7() is not new; 
however. Theorem 14. 21 and 14. 81 together yield an algorithmic/effective proof of a 
slight weakening of Theorem 14. 71 This algorithmic/effective proof is a new con- 
tribution. The material in Section^amounts to a special case of a "higher rank" 
generalisation of Kedlaya's algorithm. Section O contains a description of the 
main technique used in the "deformation algorithm" . The analysis of the loss of 
numerical precision during the application of this technique is the only original 
contribution in this section; see Theorem 15. II and the discussion following it. 
Section presents our algorithmic solution to the main computational problem. 
Specifically, we assemble together the algorithmic and theoretical techniques 
developed in Sections 01 and El to address Problem 13. 71 

Section presents an explicit family of surfaces, viz., open subsets of affine 
surfaces defined by equations of the form Z"^ = Q[X, Y) under some smoothness 
assumptions. (We note that these surfaces were previously studied for differ- 
ent reasons by the author in his expository papers |22l I23| : see also the Ph.D. 
work of Hubrechts ^3.) The algorithm described in SectionjSl together with an 
auxiliary algorithm (Section l7.3.1|l . Kedlaya's algorithm for hyperelliptic curves, 
and some propositions H7.ll 17.21 and l7.3|l allow the efficient computation of nu- 
merical approximations to the action of Frobenius on the middle-dimensional 
rigid cohomology of these open surfaces; see Theorem l7.6l for a complexity esti- 
mate. In SectionjHlwe consider smooth compactifications of these open surfaces, 
and describe how one may efficiently compute the full zeta functions of these 
compact surfaces using the main result of Section |3 see Section 18.31 for com- 
plexity estimates. We have implemented this zeta function algorithm for the 
case in which the base field is prime using the Magma programming language. 
Section [51 presents some explicit zeta functions we have computed using our 
implementation. 

The author would like to make a comment regarding the original motivation 
of this work: An interesting problem when calculating zeta functions using 
rigid cohomology is establishing good bounds on the loss of numerical precision, 
i.e., quantifying the divisions by the characteristic p which occur during the 
algorithms. It was the author's attempt to prove such precision-loss bounds by 



3 



induction on the dimension using a deep theorem of Christol-Dwork (see |S] or 
[TUl Chap. V]) which lead him to consider a recursive approach to computing 
zeta functions. The Christol-Dwork theorem, which can be thought of as a 
special case of an effective p-adic local monodromy theorem, remains an essential 
ingredient in the theoretical analysis of the algorithm presented in this paper. 

Acknowledgements: The author would like to thank his colleagues in the 
Mathematical Institute and Hertford College, Oxford, and his family for their 
support and help. He has also been greatly assisted by the kind help of Francesco 
Baldassarri, Gilles Christol, Jan Denef, Bas Edixhoven, Johan de Jong, Ralf 
Gerkmann, Kiran Kedlaya, Michael Singer, Frederik Vercauteren and Daqing 
Wan. Especial thanks to Nobuo Tsuzuki for his detailed personal communica- 
tion, and permission to include part of it as Section IXBl in this paper. 

2 Varieties and zeta functions 

Let ¥q be the field finite with q elements of characteristic p, and fix an algebraic 
closure D ¥q. For each integer s > 1, let Fg» denote the unique subfield of 
¥q of order g". Let X be a variety defined over F^. For s > 1, let |X(Fgs)| be 
the number of F^s -rational points on X. 

Definition 2.1. The zeta function of X is the formal power series 



Theorem 2.2 (Dwork). The zeta function is a rational function. Specifically, 
Z(X,T) = P{T)/Q{T) for some polynomials P{T),Q{T) G 1 + rZ[T] with 
gcd(P,Q) = l. 

Suppose that the variety X can be specified using I{X) bits of data (input 
size) and that zeta function P{T)/Q{T) requires at most 0{X) bits of data 
(bound on output size). Let S{X) = max{/(X), 0(X)} (problem size). 

The central problem in the "algorithmic theory of zeta functions" is : 

Problem 2.3. Given X compute Z{X,T) in a polynomial number of bit oper- 
ations in S{X). 

The proof of Dwork's rationality theorem can be transformed into an al- 
gorithm for Problem 12.31 see j24| . For a hypersurface the running time in bit 
operations is polynomial in (pS'(X))'^™'^'^\ i.e., it solves Problem 1 for hyper- 
surfaces assuming the dimension is fixed and the characteristic "small". This 
algorithm though is of little practical interest. The algorithm of Schoof-Pila 
solves Problem 12.31 for smooth plane projective curves of fixed degree ^51 1?7) : 
this is the most general result obtained so far using the /-adic theory. 

Let us for the remainder of this section assume that X is smooth. Then 




2dim(X) 



Z{X,T)^ n det(l-T/""WF-i|i/;,g(X))(-i) 
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The rigid cohomology spaces H^^g{X) are finite dimensional vector spaces over 
K, the unramified extension of Qp of degree e := [F^ : ¥p]. The Hnear map 
Fq is that induced by functoraUty from the qth power Frobenius endomorphism 
of the coordinate rings of an affine cover of X. We have Fq = F^ where Fp is 
semi-hnear with respect to the Frobenius endomorphism of the p-adic field K. 
We denote F := Fp. The central problem in the "algorithmic theory" of rigid 
cohomology is: 

Problem 2.4. For each i, compute a sufficienty good numerical approximation 
to a matrix for the semi-linear map F : F[^^g{X) — > H^^^^X). Moreover, do so 
in a number of bit operations polynomial in pS{X). 

By "sufhciently good" we mean good enough to recover the integer polyno- 
mials P{T) and Q{T) from their numerical (p-adic) approximations. 

In principle, Kedlaya's algorithm 18 can be applied to Problem l2.4l but the 
running time of this approach is polynomial in [pS {X))'^™^'^\ However, it is a 
remarkably useful algorithm for the case dim{X) = 1 where it has been exten- 
sively studied and implemented, see figj . Problem l2.4l was solved for smooth pro- 
jective hypersurfaces in (21) (the "deformation algorithm") using relative rigid 
cohomology. This approach again seems to be of some practical interest, see 
in particular the recent work of Gerkmann 12 and Hubrechts |13| . (We refer 
the reader to Tsuzuki 29 for a different approach which also uses relative rigid 
cohomology. This method, which one might call the "degeneration algorithm" , 
is conceptually very nice; however, it has only been worked out in one special 
case and it is not clear to the present author how widely it can be applied.) 

To understand better the performance of the "deformation algorithm" it 
is necessary to look more carefully at the dependence on input/output size. 
Specifically, one can consider separately the dependence on the arithmetic size, 
measured by \ogp{q) and p, and the geometric size. In the case of smooth 
projective hypersurfaces, the latter is the middle betti number ft,2, say, which 
is approximately {d — 1)" where d is the degree of the hypersurface and n 
is the dimension. The "deformation algorithm" has good dependence on the 
arithmetic size. However, the time/space dependence on the geometric size is 
rather high. Specifically, based on the analysis in ll^l . the author conjectures 
that the "deformation algorithm" requirements 0{plog{q)^h2^'^) bit operations 
and O{plog{q)'^h2) bits of space for a smooth hypersurface with middle betti 
number ft,2 defined over Fg. Here the Soft-Oh notation ignores logarithmic 
factors Def. 25.8], and ui is the exponent for matrix multiplication. 

The aim of the new approach in this paper is to try and reduce the space / time 
dependence on the geometric size by using a more economical geometric method. 
This is achieved for the surfaces studied in Section [T] and |S1 see Section 

We conclude this section by mentioning two very recent advances in the 
area: First, work by Kedlaya et al 1 on bounding Picard numbers using p- 
adic cohomology. Second, forthcoming work of Edixhoven and collaborators on 
computing coefficients of certain modular forms using /-adic cohomology of high 
dimensional varieties. 
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3 Rigid cohomology 



Many authors have contributed to the theoretical development of rigid cohomol- 
ogy, most notably Berthelot, Dwork and Monsky, and more recently Kedlaya 
and Tsuzuki. We follow the definitions given in Gerkmann \V2l Sec. 3] and refer 
the reader to that source for further details. 



3.1 Relative rigid cohomology 

Let k — ¥qhe SL finite field of characteristic p, and K be the unramified extension 
of Qp of degree [k : Fp]. Let Ok denote the valuation ring of K. Then (p) is the 
maximal ideal of Ok and Ok/{p) — F^- Let ordp denote the p-adic valuation on 
K normalised so that ordp(p) — 1, and | • |p := ) the corresponding norm. 

Extend the norm and valuation to polynomial rings and finite dimensional vector 
spaces over K in the obvious manner. 

Let X be a fc-scheme of finite type. Let {X,X,X) be an O^f-triple for X, 
viz. an open immersion j : X X into a proper fc-scheme, and an "admissible" 
embedding i : X X into a formal O^f -scheme. For S a fc-scheme and (S*, S, S) 
an O/f-triple for S, a morphism {X, X, X) — > (5, 5, S) is a commutative diagram 

X ^ X ^ X 

fi n if 

The relative rigid cohomology sheaf of the morphism f : X S is 

'^ligi^/S) ■■= K'fK*j^^'x[/]S[- 

We take global sections to give the relative rigid cohomology spaces H^^g{X/S) := 
T{]S[,ni.,giX/S)) with which we shall work. 



3.2 A relative comparison theorem 

Let X and S be fc-schemes of finite type. Assume now that there exist commu- 
tative diagrams 

X ^ X S ^ S 

t t and t T 

X X s ^ s. 

Here X, X, S, S are O/f-schemes. The vertical maps take the special fibres. As- 
sume that the lower horizontal maps are open immersions and their codomains 
X and S are proper and smooth O/f-schemes. Now take X and S to be the p- 
adic completions oiX and S, respectively. Then we obtain O^-triples {X, X, X) 
and {S, S, S) . Assume now that we have morphisms f : X ^ S and f : X ^ S 
as in Section 13.11 and further morphisms X S and X ^ S. Assume more- 
over that these fit together with the two diagrams immediately above to give a 
"commutative cube" . (By "completion" this yields a morphism of O^-triples 
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as before.) One may define from this commutative cube the relative de Rham 
cohomology sheaf of the induced morphism on the generic fibres : Xk —^ Sk- 

Taking global sections we define the relative de Rham cohomology spaces 
HIi^{Xk/Sk) ■= T{Sk / Sk)) with which we shall work. 

Finally, assume that the complement X — X has smooth components with 
normal crossings over S. Then TL\i^{Xk / Sk) is coherent.-^ Moreover, according 
to Gerkmann ^1 Eqn (8)], the comparison theorem of Baldassarri-Chiarelletto 
extends to this relative situation. Specifically, the natural morphism 

re^niXK/SK) ®os^ -> K^MIS) (i) 

is an isomorphism. Define := r(],5[, j^Ojgj) and A :— T{Sk,Osk) to be the 
rings of global sections. Then 0] Prop. (2.5.2) (ii)] shows that ^ implies the 
following theorem. 

Theorem 3.1. With the assumptions stated in this section, the following iso- 
morphism holds: 

W,j,{Xk/Sk) ®a At ^ K,giX/S). 

3.3 Proper and smooth base change 

We retain the definitions and assumptions from Section im and define A^ := 
T{]S[,j'^0-^s[). Assume now that the morphism X ^ S is proper and smooth. 
For each point 7 G S* in the base denote by the fibre X — > S' at 7. The 
following base change theorem will be of importance to us ^1 Thm 3.1]: 

Theorem 3.2 (Berthelot). Base change Spec(fc(7)) — > S induces an isomor- 
phism H„giX/S) (8)^1- K{j) = H„giX^). 

3.4 Pencils of varieties 

We retain the definitions and assumptions in Section 13.21 and 13.31 i.e., we have 
a morphism f : X S, along with all the auxiliary objects and properties so 
that the comparison fTheorcm 13. 1|) and base change fTheorcm I3.2|l theorems 
hold. 

Assume now that 

5 = Pi-{7i,...,7d,oo} 

where the "/i € k are distinct. Thus S* = P^. Choose 7,; G Ok so that 7^ mod p = 
7i. Take 

S = ^Ok ^ ■ • ■ ,7d,oo}- 

^This result appears to be "well-known to experts", although the author has not been able 
to find an explicit reference for it. The point is that in such a case 'H^^j^(Xk /Sk) may be 
computed via the hypercohomology of a proper morphism applied to a relative logarithmic de 
Rham complex, which is coherent. Note that in our application in Section |7| we shall prove 
finiteness directly. 
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The coordinate ring A r{SK, C'sk ) of the generic fibre Sk is the localisation 
K[r][l/r{T)] where r n^=i(r - 7i) e ^[r]. The ring of global sections 
At := T{]S[,,pO]s[) is the weak (a.k.a. dagger) completion of K[T][l/r{T)]. 
We denote this ring hy K[T][l/r{T)]^. Its elements can be written in the form 
^■g^ ai(r)r(r)' where the coefficients ai(r) g K[r] have deg(ai) < deg(r) and 
satisfy ordp(ai) — e\i\ ^ oo as |i| — > oo for some e > 0. (Note that there is 
no lower bound on the e which occur for different elements in the ring.) It is 
convenient at this stage to give a definition which we shall need later. 

Definition 3.3. We shall say that we have given effective p-adic bounds for an 
element in a = J2ieZ '^i(^)''^(^y ^ K[T, l/r(r)]t if we are given ri,5 G Q with 
77 > such that ordp(ai) > ri\i\ + S for all i G Z. 

3.5 Cohomology in the middle dimension 

We retain the definitions and assumptions from Section r3.4l i.e., we have a pencil 
f : X S and the comparison and base change theorems hold. 

Let us now focus our attention on the middle dimension. Specifically, let n be 
the relative dimension of the morphism f : X ^ S. By the comparison theorem 
and coherence of relative de Rham cohomology, we see that H^piX^ /Sk) and 
^rigi-^ I ^) ^'"6 locally free modules of finite rank over the rings A = K[T, l/r(r)] 
and At = K[r, l/r(r)]t, respectively. We shall assume that they are in fact free. 
Let us simplify our notation now by writing 

£ := HUXk/Sk), £^ H;\g{X/S), 

so by the comparison theorem (Theorem 13. 1|) we have 

£^ = £<S)A A''. 

The free modules 8 and £t come with additional structure. Specifically, deriva- 
tion with respect to T induces a connection V on the A- module £. Let us recall 
the precise definition of a connection. 

Definition 3.4. A connection "S/ on E is a map "S/ : E £ ® Vl\ such that 
V(ei + 62) = V(ei) + V(e2) and V(aei) — ei ® da + aV(ei) for all a £ A and 
ei, 62 G f . 

Here d : A ^ ^\ is the universal derivation, which in our case amounts 
to differentiation w.r.t. F. The connection induced by differentiation w.r.t. V 
is called the Gauss- Manin connection, and the pair (f,V) a \7 -module. Dif- 
ferentiation with respect to F also induces a connection Vt : ft ^ £;t (g) ^^^^ j 
here fi^^ is the module of if-linear differentials which are continuous w.r.t. the 
j>-adic norm. The comparison theorem tells us that this is just the Gauss-Manin 
connection. Specifically, assuming one can compute a matrix for V with respect 
to some basis of £ over A, then the same matrix defines the map Vt on the 
basis of £t obtained by extending scalars. This means that the Gauss-Manin 
connection Vt on £t can be computed in a purely algebraic manner. 
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The module f ■f comes with one further piece of data, namely the Frobenius 
map : £^ — !■ £t This map is induced by the functorality of the construction 
from the pth power map on the coordinate rings of an affine cover of X. It 
is additive and semi-linear. Specifically, the construction requires the choice of 
a lifting of the pth power map from A ®k to . Let us assume that we 
have chosen the obvious lifting a : A^ ^ A'^ so that a -.V ^ and a acts on 
K as the Frobenius automorphism. Then is a cr-linear map on the module 
f^. Thus with respect to a basis for £\ it can also be described by a matrix. 
Note though that this matrix has entries in the ring A\ whereas the matrix 
for V''' referred to in the previous paragraph has entries in A. Since ^ and a 
commute up to a factor pTP~^ on the weak completions of the coordinate rings 
of an affine cover of Xk, functorality of the construction yields the following 
important commutative diagram: 

F i i Finder (2) 

The data {£^,V^,F) is called a {a, '^^) -module over A^ , or alternatively an 
over convergent F -isocrystal over S. Let us recall the precise definition. 

Definition 3.5. An overconvergent F -isocrystal {£\'^\ F) on S (a.k.a. (cr, V^)- 
module over A'^ ) consists of the following data. A (locally) free A'^ -module £\ 
an injective a -linear map F : £^ ^ £K and a connection : £t £t (g) fi^^^ 
such that Diagram 0j commutes. 

We shall denote the kernel and cokernel of the map in Diagram (j^J by 
-ff^jg (S', ft ) and H^ig{S,£^), respectively. These objects are vector spaces over 
K. By commutativity, F induces a map on each of these spaces. We note 
that the space ker(Vt) = H^^g[S,£^) is certainly finite dimensional over K. (It 
embeds in the finite-dimensional space of local solutions around any non-singular 
point; an observation of Monsky.) 

We now state our final assumption on the family X S. We assume that the 
family Xk Sk comes by extension of scalars from a smooth morphism defined 
over an algebraic number field. It follows then by the "open local monodromy 
theorem" that the connection is regular, i.e., locally has only simple poles, and 
the local exponents (see Section E!T|l are rational numbers |E1 Thm. (14.3)]. 

Under this final assumption, as well as the others already in place in this 
section, it is known that H^j^g{S,£^) is also finite dimensional (follows from [3 
Corollary 2]). We will give an effective/algorithmic proof of finitcncss under 
some simplifying assumptions (follows from Theorems 14.21 and I4.8|l . (There is 
also an older, related result due to Adolphson ^ Theorem 2, Remark p. 286].) 

We now come to the main definition in the paper. 

Definition 3.6. The space E^f^^g H}.^g{S,£^) = coker(Vt : ft ^ ft^^f^i^.^. 

By a result communicated to us by Professor Nobuo Tsuzuki, when X is 
affine this space is a term in a spectral sequence for the morphism X S. In 
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fact, for X affine we have the isomorphism H'!^^^{X) = i?2Hg' ^(Iti (EJ in 
Section IXBl 

Finally we are able to state the computational problem we consider: 
Problem 3.7. Calculate a numerical approximation to a matrix for the map 

^ ■ ^2, rig ^ ^2,rig- 

We solve this problem under the assumption that we are given as input 
suitable numerical approximations to: 

• A matrix for the connection V, and this matrix has only simple poles ( even 
"modulo p ") and prepared local exponents. 

• A specialisation of the matrix for : 5^ — > £^ at a Teichmiiller point. 
We further assume that: 

• We are given as input effective p-adic estimates for the matrix for F : 

(see Definition . 

Regarding the first input, ones expects to be able to compute this matrix ef- 
ficiently in any concrete application of the method, e.g., in the case in which 
a basis of forms for £ is known it can be done using linear algebra. See for 
example our calculation in Section 17.3.11 The assumption on the matrix for 
the connection is "locally" true by the regularity and local monodromy theorem 
|14l Thm. (14.3)], since the family Xk ~^ Sk can be defined over an algebraic 
number field. Our simplifying assumption is that there is a global basis for which 
the matrix has only simple poles with prepared local exponents. 

Regarding the second input, by the base change theorem fTheorem l3.2|l the 
specialisation, say at F = 7 a Teichmiiller point, is precisely the matrix for 
the pth power Frobenius map acting on the cohomology space Hl\g{X^/¥q{^j). 
Here X^ is the fibre of the family at 7 := 7 mod p. Such a matrix can be 
computed recursively (or by Kedlaya's algorithm in the case n = 1). 

Regarding the third assumption, again one expects to be able to calculate 
such bounds in any concrete application of the method, see Section l7.3.3l 

3.6 Leray Spectral Sequence 

This section is independent of the rest of the paper. We describe the contents of 
a personal communication from Professor Nobuo Tsuzuki to the author. Note 
that the notation in this section is consistent with |30| . but varies slightly from 
that in the remainder of this paper. 

Let _ftr be a complete discrete valuation field of mixed characteristic (0,p) 
and V and k be the ring of integers and residue field of K, respectively. Let S be 
an afhne smooth scheme of dimension m over Spec(fc) and X/ S a smooth family 
with n := rel.dim(X/5') such that X is affine. Suppose there exists a smooth 
afffine lift X/S of X/S over Spec(V) with A ^ T{X,Ox) and R = T{S,Os) 
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such that ri^^y is a free i?-module. Then one can calculate the rigid cohomology 
of X/K as 

Hr^iX/K) := H^A^j, ® f2^/y) ( =: Hl,^{X/K)), 

where is the weak (a.k.a. dagger) completion of A over V and A]^ := A''iS)vK. 
Let us define a filtration Fil* of ^a/v 

Fil" lm(At^ (g)A n'^il ®R f^^/y ^ A^K ®A ^'a/v)- 
Since f^^/y is a free i?-niodule, one has 

G4ii = ^K 0^ f^^7^®flr!«^/y. 

There exists a spectral sequence [301 Thm. 3.4.1] 

Er := ®A ^ H'^+'-{A^^ ^'a/v) = m+{^{X/K), 

where the edge homomorphism is called the Gauss-Manin connection. Since 
Ef'^ = except when < q < m and < r < n, one has 

e:^-- = ... = e^^'- = h^,^-{x/k). 

Hence the top rigid cohomology group H™:^"^ {X/ K) is calculated by the Gauss- 
Manin connection: 

H^-iX/K) = cokeT{H-{A^^®An\/j,)®Rn^J^ ^ H'^{A^^®a^\,r)®r^'R,v)- 

(3) 

4 Algorithms for reduction in 

This section is independent of Section |31 although relies on it for motivation. 
We recall the definitions we shall need in an abstract manner, stripped of their 
geometric origin. 

4.1 Definitions 

Let = Fg be the finite field with q elements of characteristic p, and K the 
unramified extension of Qp of degree [k : Fp] . Denote by K an algebraic closure 
of K. Let Ok be the ring of integers of and r (F) G Ok [L] a monic polynomial 
of degree d which is squarefree modulo p. Let A := K\r, l/r(F)] and let A'^ be 
the dagger completion of A (this is defined in Section 13. 4|) . Let £ be a free 
module of finite rank m over A and define S'^ := £ ®a Let \7 : £ ^ £ ® n\ 
be a connection fPefinition I3.4|l . Fix a basis B for £ over A and represent 
elements in £ as column vectors w.r.t. this basis. Take for Q,\ the basis element 
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dT over A. Take the basis for £ (g) Q\ to be the tensor product of these two 
bases. Assume that with respect to this choice, the connection V acts as 

V = i + • ^ = ^ f ® = A"' <E> dT (4) 
di r[i ) 

where the matrix 6(r) e Mm{OK[T]) has degree in F at most d — 1. This 
assumption ensures that the matrix for V has only simple poles, including at 
infinity. This is our main simplifying assumption. Such a differential system 
is called fuchsian. Any differential systems may in principle, after a change of 
basis, be written in this form, possibly at the expense of introducing one new 
pole. See the discussion of the Riemann-Hilbert problem in Section 5.3]. 

Let : f ^ — > f ^ (g) il^i be obtained from V by extension of scalars. We are 
interested in the spaces 

KdR coker(V :£ E^'l^^g := coker(Vt : £^ . 

The notation chosen here is to remind the reader of the "geometric origin" of 
the connections we shall actually be considering. 
For R := {-/ e Qui r{-f) = 0} note that 

6(F) _ b(F) ^ 1 
r(F) r'(F)^^F-7 

where "dash" indicates differentiation w.r.t. F. Thus the residue matrix at the 
regular singular point F = 7 e i? is b{'y)/r'{'j); the set of eigenvalues of this 
matrix, denoted Ey, arc the local exponents or local monodromy eigenvalues at 
F = 7. One checks that the residue matrix at infinity is —6^-1, the negative of 
the coefficient of F''^^ in b{T). The set Eoo is defined as the set of eigenvalues 
of —bd-i- Finally, the exponent set of V w.r.t. the basis B is 

E{V,B) i^-BooU y Ey. 

76-R 

Note that this set modulo Z is independent of the basis B. 

Definition 4.1. Let p = p{\7 ,B) be the smallest positive integer larger than any 
integer in the set E{\I,B). 

Denote hy A® dTp the vector space of 1-forms spanned by the set 

I ^- ® dF : < i < rf, 1 < j < p| U {F^ rfF : < j < p - 2} . 

Denote hy S ® dTp the iiT-vector space spanned by column vectors in A™ (g) dT 
whose entries belong to the space A (g) dTp. 
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4.2 Effective finiteness of E^'^r- 

We can now state our first finiteness theorem. 

Theorem 4.2. Let the pair {£, V) be as defined in Section \jn\ and p ~ pC^ , B) 

the positive integer from Definition \4.1\ which depends upon both V and the basis 
B for £ . Then coker(V) is generated over K by the image of the space £ ® dV p. 

Proof. We shall give an algorithm for writing an element u £ £®dTm. the form 
u = V(u) + w with V £ £ and w & £ dTp. It proceeds in two stages: First, 
simultaneous reduction of the pole orders of 1-forms at the roots of r; Second, 
reduction of pole orders at infinity. 

Let [/(r) e i4r[r]™, viewed as a column vector. We shall show that for ^ > p 
we have 

U _ „fv\ W 



, , dT^V^-jj+-j®dr (5) 

for some V,W £ X[r]" with deg(y) < d and 

deg(W^) < max{max{2d - 2, deg{U)} - d, 0}. 

Moreover, we shall give a method for computing V and W. 

We claim that there exists a unique V{T) £ A' [F]™ with deg(y) < d — deg(r) 
such that 

{-lr'Im + b)V = U modr. 
Let us assume this claim is true. Define 

X := ^-^r'I„. + b)V-U ^ 
r 

Then 

deg(X) < max{max{2d - 2, deg(?7)} - d, 0}. 

Define W := —X — V'. Then deg{W) is bounded as claimed above and one 
checks by direct computation that (O holds. 

It remains to establish the uniqueness, existence and computability of V. 
For this, we must show that the determinant of the matrix (—ir'Ijn + b) is a 
unit modulo r. Now 

{-ir'I^ + b) = -r'fii,T) 



where 



f{t,r) -.^ti^n- ^^^^ 



r'(F)- 

Now r' is invertible modulo r since the latter is squarefree. We need to show 
det{f{i, F)) is a unit modulo r — YljfzR{r~l), i-e., we must show that det(/(^, 7)) ^ 
for all J E R. But det(/(i',7)) — if and only if £ is an eigenvalue of the 
matrix 6(7)/?'' (7). Since £ > p and p is larger than any integer element in the 
set U^gfli?^ the result follows. 
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Let ;7(r) G KiT]"" as before. We shall show that if 
deg(C/) - pd=:e-l> p-2 

then 

U n W 

— ®dV^\J{VT'^) + — ®dT (6) 

for some V e if™ and W G ifp]" with deg(W^) < deg(C/) - 1. 

We shall take local expansions of rational functions around the origin. Put 

u = u,.,T'-' + u.^^v'-^ + . . . , V = + (6,_ir-i + ...). 

dl 

Here hd-i G Mm{K) is the coefficient of the monomial T'^^^ in h{T). Let V G 
i<r™ be the element such that 

We note that V exists and is unique by the assumption that the integers in the 
eigenvalue set E^o of —hd-i are all less than p, and £ > p. By direct computation 
one checks that 

U W 

— (S)dr- vfyr^ = — « dr 

rP ^ ' rP 

where deg(W^) < deg(C/) — 1. This concludes the description of the algorithm. 

We note that in implementations one should represent the numerator U in 
an r(r)-adic expansion. With such a representation, in the second stage it is 
more efficient to compute 

with V G if™ as in the preceding paragraph. □ 

Theorem 4.3. The space '^^ is a finite dimensional K-vector space. 

Proof. A basis for this space can be computed using linear algebra and Theorem 
14.21 Specifically, one computes a basis for the cokernel of the J-iT-linear map 
V : ^ A™ (g) dVp. Here C A™ is the if-space of column vectors 

whose entries have poles of order at most p — 1 . □ 

We note that the dimension of i?2 d^R '^^'^ ^e calculated explicitly. 
4.3 Effective finiteness of i?2,ris 

In this section we give an effective/algorithmic proof of the finiteness of i?2'"jg 
under certain conditions. First, consider the conditions: 

(Rat.) The exponent set E{V,B) contains only rational numbers. 
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(Prep. Rat.) The exponent set E{\/,B) contains only rational numbers; more- 
over, for each E € {i5^}^gfl U {Eoo}, */ Ai, A2 G E then Ai — A2 is not a positive 
integer. 

If (Prep. Rat.) holds then we say that the local exponents are prepared. Cer- 
tainly {Prep. Rat.) (Rat.). Also, rationality of the exponents for a regular 
connection is independent of the basis chosen. 
Next, consider the following definition: 

Definition 4.4. The connection V is overconvergent if it has a basis of local 
solutions which converge on the p-adic unit disk around the generic point t with 
\t\p = 1, c.f 03 Chap. Ill, Sec. 5]. 

We now give a alternative characterisation of overconvergence. 

Theorem 4.5. The connection V is overconvergent if and only if there exists 
an element 7 £ Oj^ with ordpr(7) = such that the differential system + 
6(r)/r(r) has a basis of local solutions on the p-adic open unit disk around 

r = 7. 

Proof. That this is necessary follows by specialising the generic solution matrix 
at r = 7. For sufficiency, we observe that ^3 Chap IV Prop 5.1] allows one 
to transfer the convergence on the open unit disk around the point P = 7 to 
the same disk around the generic point. (Specifically, change variables so that 
7 = 0, and take "a" to be the generic point t with \t\p — 1.) □ 

Definition 4.6. The "dual" connection V is defined to act as 

V : A _ ^ : A"" A"' (g> dT. 
dT r(P) 

In this case the matrix b{T)/r{T) acts on the right on row vectors. 
The p-adic condition that we shall need is: 

( 0. C.) The connections V and V are overconvergent. 



Theorem 4.7 (Baldassarri-Chiarellotto). Let the pairs (f , V) and {£\V^) 
be defined as in Section \4.1\ Assume that conditions (Rat.) and (O.C.) are 
met. Then the natural morphism (f,V) — > (£^,V^) induces an isomorphism 

pl,n pl,n 
^2,dR — ^2,rig- 

Proof. This is an application of ,3, Corollary 2.6]. □ 

We shall give an effective / algorithm proof of Theorem l4. 71 under the stronger 
assumption (Prep. Rat.). More precisely, in Theorem l4.8l we give effective bounds 
on the p-adic growth of forms during the reduction algorithm in the proof of 
Theorem 14.21 under assumptions (Prep. Rat.) and (O.C). It is easy to deduce 
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surjectivity of the morphism E^'^r ~^ ^2 ^19 from Theorem 14 .81 injectivity may 
also be easily derived using the technique in the proof of j2| Theorem 2]. We 
omit the details of the proof of the isomorphism £'2 = -E'2'rig from Theorem 
14. 81 since these are not useful to us. 

We introduce notation needed for the statement of Theorem 14.81 For each 
k with 1 < k < m denote by Cfc the element of K"'^ with 1 in positiion k and 
elsewhere. For £>p,0<j<d and 1 < k < m define 

Apply the algorithm in the first stage of the proof of Theorem l4.2l to compute 

^irAi,k) ^ J2 ' , l//'^'^'^'"^ e K[Tr, deg(T//'^^'^^'"') < d 

i=p 

such that ^('^■^J.'^) - X/^v^-'Aj^^)) =: w^'^'^'^''^' G £ » dTp. Similarly, for any ^ > p 
we may apply the algorithm in the second stage of the proof of Theorem l4.2l to 
write 

where this time 

y(oo,£.fe) e^r^-i^rfr, and ^ ^ ^^(oo,£,fe)^. ^^^^ y^o^m ^ 

i=p 

We have the following effective bounds on the growth of forms during reduction 
c.f. [13 Lemma 2]? 

Theorem 4.8. Assume that conditions (Prep. Rat.) and (O.C.) hold. Then for 
w e u!(°°'^''^)} we have 

ordp(w) > - (alogp(£ - p) + /?) 

for some effective constants a, P (li Q which depend only upon the connection V 
and the basis B, i.e., are independent of the starting form u G y(oo/,fc)|^ 

We shall make the constants a,/3 completely explicit in Note l4.lTl We note 
that Theorem 14 . 81 also holds with the same constants if one applies the variant 
algorithm for reducing pole orders at infinity given at the end of the proof of 
Theorem 14.21 Since the forms ^('''^J''^) and it(°°'^''=) span E ® dV as & iiT- vector 
space, the above theorem allows one to deduce bounds on the growth of arbitrary 
forms during the reduction algorithm. 

■^For notational convenience, in the statement of the theorem and also inequaUties J7J, 151 . 
and I16i the expression logp{£ — p) occurs, when the argument could be zero; similarly logp{j) 
where j = occurs in inequality 1131 . In these cases "logp(O)" should be understood to be 
zero. 
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The proof of Theorem 14 . 81 will be reduced by a localisation argument to that 
of giving effective bounds on the p-adic convergence of the uniform part of the 
local solution matrix to a differential system at a regular singular point. Such 
bounds are provided in Lemma l4.9l whose proof in turn relies on a deep theorem 
of Christol-Dwork-Gerotto-SuUivan TXP, Chap V], and an elementary result of 
Clark (Lemma B.lOf) . 

Proof. Since 6(r) e M„i{Ok[^]) and ordp(r(r)) = 0, from the equation "w = 
u — V(u)" we see that it suffices to prove that for v e {v^'^''^'^'''\v^°°'^'''^ we 
have 

ovdpiv) > - {a\ogp{e - p) + (7) 

for some effective constants a, /? G Q which depend only upon the connection 

V and the basis B. 

We divide the proof of Q into three steps: 

• Step 1: We reduce proving bound {Tj) to proving the local bounds ®. 
Here we need that r is squarefree modulo p. 

• Step 2: We reduce proving each local bound to proving a different 
local bound (|10|) . Here we need assumption (Prep. Rat.). 

• Step 3: Bound H10|l is deduced from an effective version of a theorem of 
Christol. This step uses assumptions (Prep. Rat.) and (O.C.). 

Recall that p is defined to be the smallest integer larger than all integers in 
the exponent set E(V,B). We note that the argument we give works for any 
"p" larger than every integer in the exponent set E{y ,B). 

Step 1: First let us consider the case that u :— u'^'^'^'i^^) for some £ > p, 
< J < d and 1 < k < m. Let us simplify the notation above by removing 
the exponent ^^{r, £, j, k)" where it occurs, i.e., v := v^^'^'^'''\ w := w'^'^'^'^'^^ and 
Vi := v^^'^''''^'^ etc. Let 7 S i? be a root of r{Y) = 0. Let t-^ = F — 7 and expand 

V locally as 

= v^jt-^ + • • • + v^^pt-P + . . . , v^^^ e if(7)™- (8) 

We show now that Q holds for v — provided that 

ordp(u^,i) > -{a\ogp{l - p) + (3) for aU 7 e i? and p < i < £. (9) 

Assume Q holds. We claim that ordp(Vi) > — (a logp(^ — p) + P) for p < i < £, 
from which ^ follows immediately. This claim can be proved by descending 
induction on i in this range. Less formally, observe that v^^i is just Vi{'y)r' . 
Since the roots of r are distinct modulo p we have ordp(Vf(7)) = OTdp{vj,i). 
Since deg(V«) < d = \R\, from © we deduce the claimed bound on ordp(V^(F)). 
Now subtract V^/r^ from both sides of ((HJ and repeat the argument for « = £—1, 
and so on. 

Similarly, assuming ^ holds for the coefficients in the local expansion of 
y ■— y(co,i,k) g^-j- in^gn^ity^ ^ffQ easily deduce that {Tj) holds for v — v^°"'^'''\ 
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It remains to establish the local bound (We omit the remainder of the 
proof for V — since it is exactly the same.) 

Step 2: Fix -f £ R and simplify notation as in Step 1. Define G{t) g 
Mm(CK(7))[[i]] so that —t~^G{t) is the expansion of 6(r)/r(r) w.r.t the local 
parameter t := — T — -f. Define H := G(0), the negative of the residue 
matrix b{'~f)/r' (j). Let the local solution matrix Y{t)t^ to the differential system 
t-^ - G{t) = be defined as in ^ Chap. Ill Prop. 8.5]. (Note that we have 
chosen our signs to be consistent with JO].) The existence of such a solution 
matrix requires the assumption (Prep. Rat.). The uniform part Y{t) lies in 
Mm(iir(7))[[t]] with F(0) = /m, and the element , which is constructed on 
Uni Page 103], satisfies the equation f^{t") = f-H" , and (t^)"^ = t'" . 

Write Y = YZ^Y.t' and Y{t)-^ = YZoZ^t\ where Y,,Z, £ M„(i^(7)). 
We shall show now that it is enough to prove that 

ordp(r,), ordp(Z,) > -(ai logp(i) + /3i), z > 1 (10) 

for some explicit ai,/3i £ Q which depend only on V and B. 

Let us assume H10|l holds. Observe that we have the local factorisation 

V ^Y{t)t" o j^ot-"Y{t)-\ (11) 

Premultiplying the localised equation V(w) — u — why {Y{t)t^)~^ , using (|ll|l . 
and then integrating, we find that 

v{t) = Y{t)t" I J t-"Y{t)-\u - w)dt + c| , (12) 

for some constant c e K{'j)"\ Bound © can now be deduce by explicitly 
integrating the righthand side of (|12|l and comparing coefficients of t^^ for p < 
i < L 

Specifically, the integrand on the righthand side of (|12f) can be written 
^ t-"a,t-('+^^+\ ordp(aj) > -(«i logp(j) + /3i) for < j < ^ + 1 - p, (13) 

for some aj G if(7)™. The lower bound on ordp(aj) comes from (IIOII and 
the integrality of u{t). Note that we do not have any bounds on ordp(aj) for 
j > i + 1 — p since these terms are affected by the unknown element w(t) 
and unknown constant c. Element H13|l may be explicitly integrated "term- 
by- term". Precisely, from the defining property of the element , ones sees 
that J t-" ajt-^^+^^+Ht = {-H-{l+l) + j + l)-4-^ajt-(^+i)+^+i, plus an 
unknown constant of integration. Recall that —H is the residue matrix. Now 
for < j < £ + 1 — p we have that — £ < — (^ + 1) + ,) + 1 < — p, and since p is 
larger than any eigenvalue of —H the inverse matrix immediately above exists. 
(For j > i + 1 — p, when the inverse does not exist the coefficient aj must be 
zero.) Next, note that {—H + i)^^ for i £ Z (when it exists) commutes with ; 
this follows from the fact that H commutes with , see ^1 Page 103]. Thus 



18 



each term on the righthand side of p2|l which does not involve the constant c 
has the form 

Y,f {t"){-H -{£+!) + ] + l)-h-"ajt-'^^+^^+^+^ = Y,{-H^£ + j)-^ajt-'^+'+^ 

(14) 

for some i,j > 0. Terms on the righthand side of (|12|l which do involve c have 
the form Y,tH"c for some i > 0. Since the lefthand side v{t) is a Laurent series, 
it follows from ^0, Chap V Lemma 2.3] that cither c — 0, or c 7^ with H a 
diagonal matrix with integer eigenvalues. Since all eigenvalues of —H are less 
than p, in either case any term on the righthand side of (|12|l involving c cannot 
effect the coefficient of for p < i < £. 

From H14|) . a lower bound on the coefficient of t^^^" for < s < ^ — p on the 
righthand side of l|12() is 

min ordp (Y,{-H -e + j^^aA . (15) 

i+j=s 

We have bounds on ordp{Yi) and ordj,(aj) for « > and < j < i ~ p, viz. 1)1 Q|l 
and (|13|l . It remains to bound ordp((-ff -i + jy^) for < j < i - p. Now 
oidp{{—H — £+j)~^) > —OTdp{det{—H — £+j)) so we must find an upper bound 
for the valuation of the determinant. Denote by Ai, . . . , A™ e Q the eigenvalues 
of -H (the residue matrix). Then det(-7J - i + j) = UZii^^ " ^ + j)- Take 
the positive integer iV to be a lowest common denominator for the A^ and define 
p-i = NXi E Z; note that gcd(p, N) = 1 since the eigenvalues are p-adic integers. 
Take the positive integer A to be minimal so that |Ai| < A for all i. Then for 
0<j<i — pwe have 

ordp(A,-^ + j) = 0Tdp{p,-N{i-j)) < \ogpi\p,\+N£) < logp(7V) +logp(A + €). 

So certainly for < j < £ — p we have 

ordp(det(-ff -£ + j))<m (logp(f - p) + logp(7V) + logp(2A + 2)) , (16) 

since p < A + 1 . From (|15|l and H16() we conclude that 

a := \2a, + , /? := \2p, + m{\ogp{N) + logp(2A + 2))1 (17) 

will certainly suffice. 

Step 3: We now establish bound (|10|l . First consider Y{t). By assumptions 
(Prep. Rat.) and (O.C.), we see that the hypothesis of Lemma 14.91 are met. 
Hence we may apply the bound in Lemma |4. 91 to our differential system t-^ — 
G{t) = 0. We next note that Y{t)~^ is the uniform part of the local solution 
matrix of the "dual" differential system t-^ + G{t) = where G{t) acts on the 
right, or equivalently, the transpose of the uniform part of the local solution 
matrix of t^^ + G(i)*'^ = with G{ty acting on the left, c.f. TO', Page 193]. 
Again by assumptions (Prep. Rat.) and (O.C.) the hypothesis of Lemma 14. 91 are 
met, and we may use that bound. □ 

The next lemma is a modest generalisation of the main theorem in [101 Chap. 

V]. 
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Lemma 4.9. Let G(t) G Mm(C'x(^))[[i]] be the local expansion of a rational 
function G{T) e Mm{K{T)) around some point i = F — 7. Let S := t-^ and 
consider the differential system S — G = 0. Assume 

1. The eigenvalues of G(0) are rational numbers, and no two differ by a 
positive integer, 

2. The solution matrix to the differential system around the generic point t 
with \t\p — 1 converges p-adically on the open unit disk. 

Let Y{t) = J^iLo ^ -^m(^(7))[M] ^6 uniform part of the local solution 
m,atrix Y{t)t'-^^^^ of the differential system S—G = 0. Then there exist ai, Pi G R 
such that 

OTdp{Yi) > -(ai logp(i) + for all i > 1. 
Proof. First, change basis by a matrix H G GLra{K{'-f)[t,t^^]) so that 

G\H] ■= ^t— — h 7i ^GTi. 

^ ' at 

is such that G[-h] (0) has eigenvalues in the interval [0, 1). By ^3 Chap. V, Prop. 
4.1] the matrix Ti. may be taken to be unimodular, i.e., ordp(7i), ordp(7i^^) > 0. 
Moreover, the degree in t oi H (degree in of H~^) is the absolute value of 
the floor of the most negative eigenvalue of G{0), and the degree in of H 
(degree in t of H~^) is the floor of the most positive eigenvalue of G'(O); this is 
easily seen by viewing Ti and 7i~^ as a product of shearing transformations. 

Let Y{t) = J^iLo ^i^^ uniform part of the local solution matrix to the 

differential system 6 — G^^y Note that G[>^] € Mm(Oj^(..y))[[t]] by the unimod- 
ularity of Ti. 

The main theorem in ^1 Chap. V Section 9] assures us that there exists 
02, /32 G M such that 

ordp(f;) > -(as \ogp{i) + /^s), i > 1. 

We comment briefly on why the main theorem is applicable: In the notation of 
|1(JI Chap V], we must check conditions TZl, 712, 713' and 714. Now 711 is true 
since the matrix G^t^j contains functions which are localisations of rational func- 
tions; 7Z2 (overconvergence) follows from assumption 2. and the unimodularity 
of 7i; 712 (eigenvalues in Zp n [0, 1)) is true by assumption 1.; 714 (integrality) 
follows since G[>^] G Mm{OK('y))- 

There exists H G GL,„(i\:[t, i^^]) such that Y{t) = n-^YU, c.f. [El Page 
163 Lines 1-6]. Moreover, the degree in t of 7i is the floor of the most positive 
eigenvalue of G(0), and the degree in ofH is the absolute value of the floor 
of the most negative eigenvalue of G(0); this follows from the argument on |10[ 
Page 163 Lines 11-21]. 

bmce ordp(7^^i) > 0, to prove the lemma we need only calculate a lower 
bound on ordp(7Y). One computes this from the equation = Y^^TiY using 
bounds on the degree in t^^ oiJi and the naive upper bound on the growth of 
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the coefficients of Y{t)-'^ and Y{t), c.f. [Tni Page 191-193]. The naive upper 
bounds are given in Lemma 14.101 Specifically, one finds certainly 

ordp(H) > -/33, /33 + 41ogp(A + 1) + 2 logj,(27V)^ (18) 

where A is such that all eigenvalues A of G(0) have |A| < A and N the lowest 
common denominator for the eigenvalues. Note that we have already observed 
deg(-i ('T^"^), deg(-i (7Y) < A. Comparing coefficients in F = Tl'^YTi we see 

ordp(r,) > ordp(f,+2A) - /33 > -(^2 logp(i) + ^2 logp(2A + I) + P2 + Ps) 

which gives a bound of the required form. Precisely, take 

ai a2, Pi aa logp(2A + 1) + /Jz + /Sa- (19) 

□ 

The following lemma is an effective version of the general bound of Clark 

ElEHI. 

Lemma 4.10. Situation as in the statement of Lemma \4-y\ only without as- 
sumption 2. Denote by Ai,...,Am the eigenvalues of G{0) and assume that 
I Ail < A and NXi G Z for all i and some minimal integers A > 0, A'^ > 1 with 
gcd(p, N) = 1 . Then for all s > 1 we have 

ordp(n) > -m' (^-^ + logp(l + s) + logp(2A + 1) + log^iN)^ . 

Proof. The power series Y{t) may be computed using the classical method in 
|l()l Chap V, Remark 2.2]. It follows from this recursive method that 

s m 

ordp(Kj) > - X! X! o^dp(fc - Sij) 

k—l — l 

where 5ij :— A; — A^. We shall estimate the righthand side using the method 
of Clark. Specifically, for a & Zp and s a positive integer, define 9{a, s) := 
Yll^iict + k) . Then ordp(Kj) > —^ijOTdp9{—6ij,s). We compute an upper 
bound for the 9{~Sij,s) using the argument in [HJ Page 265, Case 3]. First, 
write —5ij — Vij/N, where < 2iVA. Then for any positive integer s we 
have 

ordp(-,5^j+s) = orAp(v^j+Ns) < logp{\v,j \+Ns) < logp(iV)+logp(2A+l)+logp(s) 



This shows that we may take the expression "w(a;) — — fclogp(l + x) ~ k 
immediately preceding Eqn (13)], to have coefficients "fc — 1" and "fc' = 
logp(2A + 1) + logp(Af)" (we have changed Clark's "log" to "logp"). From 
Eqn (14)] we deduce that 

ordp0(-%, s) < ~ vis) = + logp(l + s) + logp(2A + 1) + logp(A^) 

as required. □ 
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Note 4.11 The constants a,(3 £ Q in Theorem l4.8l can be made completely ex- 
plicit. Precisely, by equations H17() . H18() and (|19l) one sees that it suffices to make 
the constants a2 and (32 from the proof of Lemma explicit . These constants 
are those which occur in the theorem of Christol-Dwork-Gerotto-SuUivan. 

The theorem of CDGS states that ordp(Ki) > -(a2 Llogp(i)J +132) with 
012,(32 € M as follows. Define 



When all the eigenvalues are zero (nilpotent case) we can directly apply the 
Christol-Dwork theorem ^1 Chap. V, Thm 2.1]. This gives a2 = Bm,p and 
(^2 = 0. In particular, forp > m we have a2 = m— 1, and in general a2 < 3m— 3. 
For eigenvalues in the interval [0, 1), one applies the generalisation in 10. Sec. 
9, Chap. V]. Define t := [logp(i)J + 1. From the equation on middle of page 
198 in JO! one deduces that 



where the number B is defined on |10l Pages 197-198]. (Note that we have 
changed from multiplicative to additive notation, so our "_B" corresponds to 
"logp(B)".) Let N be the lowest common denominator for the eigenvalues of 
the residue matrix. Then for p > {m, 2N} from ^Oj Remark 2.2] one sees that 
Bm.p — m — 1 and _B = 0, so a2 = — 1 and P2 = 2m^ — m — 1. For general 
p, one computes from jlOJ Page 198, Line 7] and ^1 Page 197, Line 12] that 
B < {£+ l)m{m — 1) logp(2A'^), and as observed before Bm,p < 3m — 3. Thus 
we may take 

as = m(m-l)(l + logp(2iV)) + 3m-3, (32 = 2m(m - 1)(1 + logp(27V)) -f 3m- 3. 

We conclude the following: if A > is a bound on the absolute value of the 
local monodromy eigenvalues, > 1 a lowest common denominator, and m the 
dimension, then one has 



for some absolute constants Ci,C2 G Q- This will be useful in our complexity 
estimates. 

In Section [6.21 we shall see that conditions (Rat.) and (O.C.) are met in 
the situations (which arise "from geometry") which we shall encounter. The 
stronger condition (Prep. Rat.) will be met in the examples we compute in 
Section El 

Theorem 14. 81 is essential in the complexity analysis and practical application 
of the algorithm in the proof of Theorem 14.21 for the following reason: When 
one calculates the reduction of differential forms using this algorithm it is im- 
practical to store the coefficients "exactly" . At each step of the reduction, one 




(20) 



ordpiYi) > - {(^ + l)(m(m - 1)) + W,n.p + B} 
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"approximates" the coefficients modulo some fixed power of tlie characteristic. 
Making this approximation amounts to adding an "error form" to the form 
being reduced. Theorem 14.81 shows that the error introduced propagates in a 
"logarithmic" manner during the remainder of the computation. Furthermore, 
Theorem 14 . 81 applied with a general value "p", at least as big as p itself, allows 
one to bound the intermediate coefficient size during the reduction computation. 
Note that a naive inspection of the reduction formulae in the proof of Theorem 
14.21 suggests that terms grow and errors propagate in a "linear" manner; that 
this is not the case for calculations in rigid cohomology was an important insight 
of Kedlaya c.f. ^1 Lemma 2]. 

A similar "logarithmic error propogation" phenomenon arises in the numer- 
ical solution of differential systems, as we shall see in Theorem 15. II of the next 
section. 



5 Deformation of Frobenius 

In this section we retain the notation and definitions in the first paragraph of 
Section Wa\ but slightly alter some of our assumptions. Specifically, f is a free 
A-module of rank m, where A = K\r^ l/r(r)] and K is the unramified extension 
of Qp of degree [F^ : Fp] . The map V : £ ^ S ® ^\ is a connection which with 
respect to a fixed basis of £ (and "natural" corresponding basis oi £ ® Vt\) acts 
as: 

Here b{T) £ A/„i (O/f [F] ) . In this section we do not assume that the connection 
has only simple poles, i.e., we do not need the assumption that r(F) e Oifp] is 
squarefree, nor do we need the degree restriction on the matrix b{T). However, 
we shall add the new assumption that r(0) ^ mod p. 

Let : £'^ £'^ ® fi^^ be obtained from V by extension of scalars; the ring 
At = K\r, l/r(F)]t is described explicitly in Section El Let a : A'^ he 

the lifting of the pth power Frobenius which maps F ^ . Let F : £'^ £'^ he 
an injective a-linear map such that the triple (ft, V^, F) defines a (ct, V)-module 
over A'^ (see Section IS75|) . i.e., the diagram ^ commutes. 

5.1 Local deformation 

If we assume that the connection matrix B{T) is known and that the specialisa- 
tion F(0) is also known, the commutative diagram (0) allows the computation 
of a local expansion of the Frobenius matrix F{T) around the origin to any 
required precision. We describe two different approaches. 
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5.1.1 Method 1 



Let C(r) be a basis of local solutions to the differential system V = with 
initial condition C(0) = Im- So 

rlC 

— + B{r)C{T)^0. (22) 

Commutativity of l(2Jl implies that the Frobenius map F is stable on the kernel 
of the connection. Recalling that the map F is cr-linear, we deduce the matrix 
equation 

F(r)C""(rp) = c(T)D 

where D is some constant matrix. Evaluating both sides at F = shows D = 
F{0). So we have the local factorisation 

F(F) = C(F)i^(0)(C"^(FP))"i. (23) 

Thus F(r) can be computed modulo F^"" for any Nr > 1 provided we can 
compute C(F) modulo F^"". A simple recursion formula for computing the 
matrix coefficients in the local expansion of C(F) — J27Lo ^i^^ '^^^ be derived 
from the equation 

rlC 

r(F)— +6(F)C7(F) = 0. 

Specifically, write r(F) = J^Zo"^^ ^^r^ and 6(F) = E£o ^<r\ Then for ^ > 1 
we have 

/ dcg(f)) dcg(r) \ 

Ci^-—Aj2 ^^c-d-i)-, + ^'(^ - *)^^-^ j ■ (24) 

One can, of course, compute the series C°^(F^)^^ by power series inversion. 
However, it is better to observe that the matrix C(F)^^ is the solution of the 
"dual equation" 

- C{T)-'B{r) = 0, C(0)-i = J™. (25) 

It is impractical to carry out the above computations using "exact arith- 
metic"; one desires to "truncate" each coefficient Cg "modulo p^" for some 
appropriate > after it has been computed. It is an essential task to analyse 
the "propagation" of the error this introduces as one continues the computation. 

Let Eg e M,n{OK) for £ > 1, and be a non- negative integer. Let the 
sequence Di e M,n{K) for £ > be computed in the following manner. Define 
Do ■= Im and for i> I, 



^ / dcg(f)) dcg(r) 

\ i—0 2—1 



l+p^^;,. (26) 
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The series D{T) := X^fco ^^^^ an "approximate solution" to the differen- 
tial system (|^ computed "modulo p^" . In practice, the "error sequence" Ei 
is chosen to ensure the p-adic expansions of the entries of De are "truncated 
modulo p^" . 

Theorem 5.1. Let C{T) be the solution to ^ with C(0) = /,„, and D(T) the 
"approximate solution" defined via equation l^2b\l . Then for i > we have 

ovdpiDi -Ce)>N^ a'Llogp(^ + 1)J 

for some explicitly computable constant a' > 0. Furthermore, one can take 
a' = 6m — 5 for any prime p, and a' = 2m — 1 when p > m. 

Proof. First observe that we have the local factorisation around the origin 

V = Cir)o±oCiTr\ (27) 

Next observe that the series -D(r) is a local solution to the inhomogeneous 
differential equation 

fin 

r-(r)^ + b{r)D{T) = p^ro J] iEeT'. 

1=1 

Thus y{D) = p'^ E{T)r{T)~^ where E{T) := roJ2T=i^EeT^- Using the local 
factorisation H27|l one deduces that 

^ {cirr^DiT)) = c(r)-V"i?(r)r(r)-i. 

Integrating we find that there exists a constant matrix c such that 

D(r) = c{T) (^J c{r)-^p^ E{T)r{r)-^dr + c^ 

Note that E{T)r{r)-^ e rOA'iir]] since r(0) is a unit. Since C{T) = D{T) mod 
r we deduce that c = Im- Hence 

D{V)-C{T)^p^C{V) J CiT)-^EiT)r{T)-^dT. 

The connection V and its dual V come from overconvergent f -isocrystals, viz, 
(f^,Vt,F) and its "dual V^^, F~^)" . Hence Dwork's trick of analytic con- 
tinuation via Frobenius |15[ Prop. 3.1.2] shows that condition (O.C.) is met. 
Moreover, the connections are regular at zero, so local exponents are all zero. 
Thus we can apply the Christol-Dwork theorem to deduce effective logarithmic 
bounds on the growth of coefficients of C(r) and C(r)~^. Moreover, integration 
only has a "logarithmic" effect on the growth of coefficients of a power series. 
Explicitly, we can use the constant in the original Christol-Dwork theo- 

rem, see (|20ll in Note 14 .111 to deduce a' — 2B,n,p + 1 and the constant "/?' = 0". 
This completes the proof. □ 
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Let Np and N be positive integers. Let -D(r) be an "approximate so- 
lution" computed "modulo p^" to the differential system (|22|l modulo F^^. 
Let D{T) be an "approximate solution" computed "modulo p^" to the dual 
system ^ modulo T^^^/p\ Let G(0) G A/™ (if) be such that ordp(i^(0) - 
G(0)) > p^. Define G(r) := D{r)G{0)D'' (TP) modF^r. This is our approx- 
imation of the local Frobenius matrix F(r). We need to bound from below 
ordp {{F(T) mod F^r) _ G(F)). 

From Theorem lO L>(F) = G(F) -f p^'e(F) mod F^^ where ordp(e(F)) > 
withAf' := N-{2B^^p+l)l\ogp{Nr)\, and D"^ (TP) = {TP))-^+p^" e{r) mod 
F^r where ordp(e(Fj) > with N" :^ N ~ (2B„,p + l)llogp{\Nr/p])\. Note 
that from the Christol-Dwork theorem we have 

ordp (G(F) mod F^r) > -B^^pl\ogp{Nr - l)\ 
ordp ((C-(FP))-i mod F^^) > ~B,n.pl^ogpi\Nr /p] - 1)J. 

One now readily calculates a lower bound on 

ordp (^{CiT)F{0){C''iTP))-^ mod F^^) - L'(F)G(0)L''"(Ff )) 

to be 

min {N' + ordp(F(0)) -|- ordpiC (TP)-^ mod F^^-), 
N + ordp(G(F) mod F^^) + ordpiC (TP)-^ mod F^^), 
N" + ordp(G(F) mod F^^) + ordp(F(0))} 

>N- (3B„,,p + l)Llogp(A^r)J + B„,,p + min{ordp(i^(0)), 0} . (28) 

For example, when F{0) has integral entries and p > m we have that the loss of 
accuracy when computing F(F) mod F^"" is bounded by (3m — 2) [logp(Afr)J — 
(m- 1). 

5.1.2 Method 2 

The approach in this section is based upon that taken by Tsuzuki |29j. We do 
not give an analysis of the propagation of errors for this method, although this 
is an interesting problem. 

Commutativity of diagram ^ implies that 

HF 

— + B(F)i^(F) = pTP-^FinB^iTP). 
dT 

For £ > 1, the coefficients Fi in the local expansion F(F) — X^fco ^i^^ '^^^ be 
found recursively by rewriting this equation in the form 

dF 

r(F)r"(FP)— +r"(F?')6(F)F(F) =pFP-V(F)F(F)6"(F'') 

and equating the coefficient of F^^^ on both sides. This more direct method 
eliminates the multiplication of power series needed to compute the righthand 
side in (|23|l and is also more space efficient. 
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5.2 Global deformation: analytic continuation 

The entries in the Frobenius matrix i^(r) are p-adic holomorphic functions on 
the p-adic projective hne with open unit disks around the poles of r removed, i.e., 
uniform hmits of rational functions on this closed domain Di , say. Therefore, 
they can be uniformly approximated on this domain Di modulo any power of 
p by a matrix of rational functions whose denominators are powers of r{T). 
Using the method in Section one can compute the local expansions of these 
holomorphic functions to any required p-adic and F-adic accuracy. We now 
sketch how to "analytically continue" these local expansions, i.e., how given 
a power of p one can compute the rational functions which approximate the 
entries in the Frobenius matrix to that power. 

The essential point is that the theory guarantees that the holomorphic func- 
tions in the Frobenius matrix F(r) are "overconvergent" . This implies that they 
converge on the p-adic projective line with open disks of some unknown radius 
s < 1 removed around the poles of r. Let us notate this unknown larger domain 
by Ds- Assuming one has an upper bound on s, and also an upper bound on the 
maximum value t taken by the Frobenius matrix on the closed set Dg, one can 
compute an upper bound on the total degree of the rational functions needed 
to approximate F{T) on this domain modulo any given power of p. This upper 
bound allows one to determine how many terms in the local expansion of F{T) 
are required to compute the rational functions. The knowledge of bounds s and 
t amounts to having effective lower bounds on the p-adic decay of the entries 
in the matrix F(r) fPefinition 13.3(1 . We shall assume that these effective lower 
bounds are known; for the explicit example we consider in Section [3 we will 
explain exactly how to calculate them. 

We refer the reader to |21[ Section 8.1] for a detailed description of the rel- 
atively straightforward step of recovering the matrix of approximating rational 
functions from the local expansions given that these bounds are known. 

6 An algorithm for computing F : E^'^^^ ^I'rig 

In this section we gather together the results from Sections |21 El and [S] and 

present the main algorithm of the paper. 

6.1 Definitions and assumptions 

In this section we retain the definitions from Section [3. 51 and make the assump- 
tion on the connection matrix from Section 14.11 Specifically, we are given a 
pencil X ^ S oi fc- varieties such that: 

• The relative space ft Hl\g{X/S) H2r{XkISk) ®a is a free At- 
module of rank m. 

• The connection : ft ^ (g, j^i^^ jg given by a matrix 6(r)/r(F) with 
simple poles; here r{T) € Ck[F] is squarefree modulo p and d :— deg(r). 
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• The space Hg ^ '^^ Definition 13.61 

• The morphism X S is proper and smooth, so that the base change 
theorem holds fTheorem I3.2|l . 

• The morphism Xk Sk arises by extension of scalars from one defined 
over an algebraic number field, so the local exponents are rational. 

6.2 The comparison theorem in the geometric case 

Theorem 6.1. With definitions and assumptions as in Section \6.1\ Condition 
( 0. C.) is met. 

Proof. We have an overconvergent F-isocrystal (f'',V^,i^) on . Choose any 
point r = 7 G Ojf such that ordp(r(7)) = 0. Dwork's trick of analytic continu- 
ation via Frobenius |15l Prop. 3.1.2] tells us that the basis of location solutions 
to the differential system V = converge on the open unit disk around F = 7. 
The same is true for the second differential system V = 0; consider the "dual 
F-isocrystal , V'^ , F~^)" and apply Dwork's trick once again. So by Theorem 
14. 51 condition (O.C.) is met. □ 

Theorem 6.2. With definitions and assumptions as in Section ld. 11 condition 
(Rat.) is met. 

Proof. Since Xk Sk can be defined over the complex numbers, it follows 
from the local monodromy theorem 0%!, Thm (14.3)]. □ 

Theorems 14.71 16.11 and 16.21 together yield: 

Theorem 6.3. With definitions and assumptions as in Section \6.1[ . the com- 
parison theorem E^ '^r — ^2 rig holds. 

We note that our bounds on the growth of forms fTheorem l4.8|l only holds 
when the assumption (Rat.) is replaced by the stronger assumption (Prep. Rat.). 
This will not always hold in the geometric case; however, we will give examples 
in which it does hold (Section O. 

6.3 Numerical approximations 

In this section we formalise the notion of a "numerical approximation" to a 
p-adic number. 

We assume that elements in Ok are represented as p-adic expansions with 
coefficients in some fixed set of representatives for the quotient Ok/{p)- Thus 
for any positive integer N, elements in the quotient Ok/{p^) can be rep- 
resented in a unique manner via truncated p-adic expansions. Elements in 
Ok[^, l/^(r)]/(p^) have an obvious representation via these truncated p-adic 
expansions. 
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Definition 6.4. Let N be a positive integer, and a & = K[r, l/r(r)]'l'. De- 
fine N' A^— inin(0, ordp(a)). A -approximation to a is a triple {N, ordp(a), ao) 
whereao G Oi^[r, l/r(r)]/(p^') andao-p-"'"(°'"'^''(°)'°)a = mO^ir, l/r(r)]/(p^'). 

Thus taking So to be any preimage of oq in Ok[^, ^/r{T)] (g) K, we have 
that a p^-approximation to a defines an element ai :— p™'"(°''^p(°)'0)ao such 
that ordp(a — oi) > iV. Conversely, given such an element ai one may canoni- 
cally identify with it a p-'^-approximation to a. Intuitively a p^-approximation 
amounts to knowledge "modulo p^" . 

6.4 Input/Output specification for the algorithm 

We retain the definitions and assumptions from Section Ki.ll and now further 
assume that (Prep. Rat.) holds. So by Theorem 16.31 our comparision theorem 
^I'dR — ^I'rig holds, and moreover, we have effective bounds on the growth of 
forms during reduction f Theorem 14. 8|l . 

Let Nj be a positive integer. We shall assume that we are given as input the 
following: 

• Input 1: The matrix 6(r)/r(r) for the connection V. 

• Input 2: A -approximation to -^(7) for one Teichmiiller specialisation 
r = 7 of the matrix F(r) for the action F : Hl\g{X/S) Hl\g{X/S), 
i.e., an approximation to the pih power Frobenius action on Hl\g{X^) for 
some fibre of the family X S. 

We also assume we are given (see Definition 13 .bfl : 

• Input 3: Effective p-adic bounds on the entries in the matrix F(T). 
The algorithm gives as output: 

• Output: A -approximation to a matrix for F : -^2 ^9 ^ ^I'rig^ 

for some effectively computable Nq < Nj. The loss of accuracy is measured by 
the difference Nj — Nq- We note in Section IHTI that there exist effectively com- 
putable constants a" , (3" > such that one may take Nq := Nj—{a" \og^{Nj) — 
/5"), i.e., we have a "logarithmic" loss of accuracy. 

When X is affine we have from equation ^ that iJ"j^^(A") = E^'^ig and 
so the matrix given as output yields an approximation to pth power Frobenius 
action on (AT). 

6.5 The algorithm 

The algorithm comprises two steps 

• Step 1: From Inputs 1,2 and 3 use the "deformation algorithm" to compute 
a p^-approximation to the matrix for F : H^\g{X/S) H^^g{X/S). 
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• Step 2: From the output of Step 1 and Input 1, use the algorithm from 
Section2]to compute a -approximation matrix for F : i?2'rig ~^ ^I'rig- 

The intermediate precision N, where Nq < N < Nj, can be computed from the 
input data, see Section [6.7l 

6.5.1 Step 1 

This step was described in detail in Sections l5.1l and l5^ (One may need to make 
a change of basis so that 7 = 0.) Estimates from Input 3 determine the F-adic 
accuracy required in the local deformation to compute an p^-approximation to 
the global matrix F(r) itself, whose entries lie in . 

6.5.2 Step 2 

Let ;8 be a set of forms in A'"^ (g) dT whose image in £'2'"ig ^ if-basis for this 
space. This may be determined by computing the set of exponents i?(V, B), then 
performing the linear algebra computation described in the proof of Theorem 
14.31 Theorem l4 . 71 assures us that this set gives a basis for El'r^g. We can assume 
that ordp(e) = for all e £ B. 

For each e G B, one computes a p^-approximation of the image 
F{T)pTP^^e'^{TP)dT. Then one performs a radix conversion from F-adic to 
r(F)-adic expansions, so that the input is in the appropriate form for the re- 
duction algorithm. (It is actually much better in practice to compute an r-adic 
expansion of the matrix F (T) / r'^ (TP) , and recover the r-adic expansions for each 
image form via a r-adic multiplication routine. It turns out that these radix 
conversions are very time consuming, so one wishes to minimize the number 
performed.) Then use the reduction algorithm from the proof of Theorem 14.21 
plus the final linear algebra step from the proof of Theorem 14.31 to write this 
as a -approximation to a if-linear combination of elements in B plus an 
-approximation to an element in V^{£^). 

6.6 Analysis 

6.7 Loss of accuracy 

Theorem 14.81 and the analysis following Theorem 15.11 together show that the 
matrix computed in the algorithm is a p^'' -approximation to a matrix for the 
action of F : i?2,'"ig ~^ ^'^'^ such that 

No^Ni-{a"\ogp{Ni)+l3"). 

Here a", (3" > are constants which may be computed from m, p, the local 
exponents of the connection, and the effective p-adic bounds on F{T) (Input 
3). We shall not present explicit formulae for a" and (3" in the general case as 
they are rather complicated. Wc note that the discussion following Theorem 
15. II allows one to compute the intermediate precision which is attained after 
Step 1. 
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6.8 Time and space complexity 



The time and space complexity may be calculate given the effective p-adic 
bounds on the matrix F{T), and also a bound on the height of the local mon- 
odromy eigenvalues, c.f. Section [7. 5. 21 We do not present an explicit expression 
for the general case since it is rather complicated. Let us just make a few ob- 
servations on Step 1: The calculation of the local solution matrix in Step 1 
(Sect ion l??T|l is fast, both in theory and practice, since it just requires the itera- 
tion of a short linear recurrence; however. Method 1 is rather space consuming in 
comparison to Method 2. The analytic continuation step requires only a single 
multiplication by a power of r(T) (computed modulo a power of F). The radix 
conversion, though in theory "quasi-linear time" 11, Alg. 9.14], is in practice 
rather time consuming. 

7 The Frobenius matrix of an affine surface 

In this section we apply the algorithm in Section to compute to any required 
numerical precision a matrix for the pth power Frobenius map acting on the 
middle dimensional rigid cohomology of a certain affine surface. Specifically, we 
consider an open subset X of the affine surface defined by an equation of the 
form = Q{X,T), subject to certain smoothness assumptions. The algorithm 
from Section 16.51 allows the efficient computation of an approximation to the 
Frobenius map F : H^^g{X) H^^g{X), provided one can obtain the auxiliary 
inputs 1, 2 and 3 (Section 16.4(1 . After defining the surface in Section ITTl we 
describe how the necessary auxiliary inputs may be calculated (Section 17. 3|) . 
Having specified some local monodromy restrictions to ensure applicability of 
the algorithm in Section (see Section I7.4|l , we then give a precise complexity 
analysis (Theorem 17. 6|) . 

In Section IHl we shall apply the results of the present section to compute the 
full zeta function of a compactification X of the open surface X. We report on 
our Magma implementation of this final algorithm in Section IHl 

We retain the notation in Section |31 In particular, recall that A; = is the 
finite field with q elements, and K the unramified extension of of degree 
[k : ¥p]. We assume now that the characteristic p is odd. The ring of integers 
of K is denoted Ok- Let us further assume we are given L D Q an algebraic 
number field with ring of integers O^. Assume that we have an embedding 
Ol C Ok and Oi/(p) = ©^/(p) = F,. 

7.1 Definition of the pencil 

Let QiX,T) e OL[X,r] and denote by Q(X,F) e k[X,T] its reduction modulo 
p. We shall assume that both Q and Q are monic in X of degree 2g + 1 where 
gcd(p, 2g+l) = 1. Let f (F) := Res(X, Q, |§) G Cl[F] be the Sylvester resultant 
w.r.t. X of Q and see Pages 150-151] or Section ITXTl (The notation 
r(F) is reserved for the monic denominator of the connection matrix, which is 
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a factor of r, see Section [^O]) Assume that f (F) has leading coefficient a 
unit modulo p, and f(0) ^ modp; in particular, it does not vanish identically 
modulo p. Define the Oif-schemes 

X Spec(B) where B := Ok[X, T, Z, f(r)-'^]/{Z^ - Q{X, T)) 

and 

S := Spec(A) where A := ©^[r, l/f(r)]. 

Let B := B ®Ok ^ ^-^d A := A ^ be the reduction of the coordinate rings 
modulo p. Define the fc-schemes 

X := Spec(i?) and 5 := Spec(A). 

We have the obvious commutative diagrams 

A ^ B X ^ S 

I i and i i 

A B X ^ S 

where the vertical maps in the second diagram take special fibres. Recall that 
the generic fibres are denoted Xk and Sk, respectively. The horizontal maps in 
the second diagram are smooth morphisms of smooth schemes, and the fibres 
are (affine) hyperelliptic curves. 

The relative cohomology spaces which concern us are: 

£ := H}j,{Xk/Sk) = ( \0<i<2g 



At 



Q 

and 

£t ■.= Hy^{X/S) = (-J=-\Q<^<2g 

Here \/Q denotes the image oi Z in B (precisely, in B^ for the second space) . We 
refer the reader to Section rOI for a description of the ring A"^ . That H\j^{X /S) 
and Hl^g{X/S) are spanned by these forms is shown in [22 Sees. 4.2, 5.4]. That 
they form a basis follows by a specialisation argument, and the fact that the 
dimension of the first de Rham (rigid resp.) cohomology space of any fibre in 
the family X ^ S {X ^ S resp.) is 2g. Alternatively, see ^S]- Note that we 
do not need to appeal to the finiteness and comparison theorems in Section [3. 21 
since we can establish the necessary results directly. 



7.2 The spectral sequence 

The next proposition shows that the algorithm in Section in the present case 
computes a numerical approximation to a matrix for F : H^.^g{X) H^^g{X). 

Proposition 7.1. With the morphism X S as defined in Section \ 7.1\ we 
have 

Proof. Follows from equation Q since X is affine and ri^^y is free of rank one; 
here R F(5, Os) and V := Ok- □ 
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7.3 Auxiliary data: Inputs 1,2 and 3 



We now explain how to compute the auxihary information needed as input 
to our main algorithm (Section |SJ) in the case of the surfaces presently under 
consideration. 

7.3.1 Input 1: the matrix for V 

The connection V acts on £ by taking the derivative w.r.t. F of the basis 
elements X'^dXj \fQ^ and then applying the reduction algorithm of Kedlaya 
to write the image as a linear combination over A of the basis elements plus 
■jx{,g)®<S^ for some element g & £ ■ We now give explicit formulae for computing 
the matrix for the connection (based upon Magma code written by the author). 

For an element a e iy(r)[X] and « e Z, denote by Coeff(a,z) the coefficient 
of X'^ in a. Let (5 := 2g + 1. Let M be the Sylvester matrix w.r.t. X of Q and 
|§. Explicitly, M e M2S-iiL[T]) and for 1 < j < 25 - 1 

_ r GoeS{X^-^-'Q, 26-1- j) for 1 < i < 5 - 1 
\ Coeff(X2*-i-'|§, 25 - 1 - i) for S<i<2S-l. 

We have assumed the determinant of this matrix f(F) (Sylvester resultant) is 
non-zero modulo p. Define E to be the S — 1 x 26 — 1 matrix over L[r] with 

E,, := -icocff I^X'-'^, {25 -I)- . 

Let F := EM~^, a5-lx 2(5-1 matrix over L[F, l/r(r)]. Let the vectors 
a, 6, c e L[F, l/f[T)][XY-^ be defined as follows: For 1 < i < 5 - 1 

ITI2 • — di ^ 2 . 

Then the connection matrix B{r) e Ms-i{L[T,l/f{r)]) is defined by Bij := 
CoeS{mj,i — 1). One can uniquely write -B(F) = 6(F)/r(r) where r\r, and r is 
monic and coprime to some entry in the matrix b G Ms-i{L[T]). 

We shall impose some restrictions on the connection matrix B{T) in Section 

7.3.2 Input 2: The Frobenius matrix of a fibre 

We take the fibre at 7 = 0, noting f(0) 7^ mod p. The Frobenius matrix of the 
fibre Z"^ = Q{X,0) can be computed using Kedlaya's original algorithm [TH) : 
the implementation by Michael Harrison is available with the documentation 
accompanying the Magma program. 
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7.3.3 Input 3: Effective p-adic bounds for F : S'^ ^ 

The Frobenius matrix F(T) can in principle be calculated by applying Kedlaya's 
algorithm to the "generic" hyperelliptic curve in the family, which is defined over 
the function field Fq(r). From the point of view of complexity theory this is 
not a good idea; it is faster to use the indirect method of the "deformation 
algorithm". However, this direct method is a good way to calculate effective 
p-adic bounds for the matrix F{T). 

Specifically, fix z, j with 1 < i,j < 2g. Let /(F) be the {i,j)th entry in the 
matrix F(r). Then /(F) is the coefficient of X^dX/ y/Q in the expression one 
obtains by reducing the form 



/XidX\ _ pXP(^+i)-i / QP - Q" [XP ,TP) 



-1/2 

dX 



using the "generic" version of Kedlaya's algorithm. Here a is the map sending 
F I— > FP, X I— > XP , and acting like the pth power Frobenius automorphism on 
K. We can write Qp - Q''{XP,TP) = pR{X,T) for some unique R e Ok[X,T] 
with 

deg^(i?) < pdegxiQ), degr(i?) < pdegr(O) 
with second inequality strict if Q is monic in F. Then 

aix^dx/VQ) = E {'Ty^'^^^s^'^- 

The ^th term in this series can be reduced modulo exact forms using pi + [2J 
applications of Kedlaya's "pole reduction formula", see |23 Section 4.2]. Each 
application requires one division by the resultant f{T). By an easy specialisation 
argument, ^18. Lemma 2] implies that reduction of the £th term requires a 
cumulative division by at most pLiogp(p(2^+i))J ^ 
Write 

oc 

/(r)= E M^^y 

k—~oo 

where fk G K[T] with deg(/fe) < deg(f). The argument in the preceding para- 
graph implies the following. 

Proposition 7.2. For k < we have the lower bound 

ordp(/fc(F)) >(£ + !)- [\ogp{p{2e + 1))J 

where £ is the smallest integer such that p£ + [|J > \k\. (Explicitly, i := 
L^^2^J and so ordp(/,(F)) > l\k\/p\ - Llogp(2|fc| + 1)J 

A lower bound for /c > requires a more detailed analysis: Let Adj(M) 
be the adjoint of the Sylvester matrix. Each application of Kedlaya's pole 
reduction formula increases the degree in F (degree of numerator minus degree 
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in denominator) by degp ( Adj (M)) — degp (r ) . The degree in F of the numerator 
XP(J+i)-ii?^ of the hh term in the series is f degr(i?) < £ipdegj.{Q) - 1). Thus 
after pi + apphcations of Kedlaya's pole reduction formula, the degree in T 
of the reduction of the £th term in the series is at most 

k{£) (pe + ^) (degr(Adj(Af)) - degr(f)) + £degr(i?). 



We note that the modest use of Kedlaya's formula for reducing the "pole at 
infinity" required in the calculation of F{r) does not increase the degree in F 
(or introduce powers of f(F) on the denominator). Thus we deduce: 

Proposition 7.3. For k > we have the lower bound 

ord,(/fc(F)) > + 1) - [\ogpipi2i + 1))J 

where £ is the smallest integer such that > k. (If no such £ exists then 

the term /a;(F) is zero.) 

Explicitly, define 

^ ^ degr(Adj(M)) ^ ^, ^ degAR) ^ degAQ) _ ,39) 
degp(f) ' ■ pdegr(f) ~ degr(r) ' 

Then assuming S + S' > 1 one takes £ the floor of ^^2p(a+a'-i)^^ • 
We now state a conjecture to which we shall refer later. 

Conjecture 7.4. The Frobenius matrix F(T) has a pole of finite order at in- 
finity, rather than an essential singularity. 

Conjecture 17.41 thus claims that /fc(F) = for sufficiently large positive k. 



7.4 Local monodromy assumptions 

In this section we state some further restrictions made to ensure that the con- 
ditions required for the application of the main algorithm in Section |S1 are met. 
Specifically, we need that the connection matrix B{r) from Section [V. 3. II is of 
the form required in Section f5.1l and that condition (Prep. Rat.) is met. To sim- 
plify the complexity analysis and to keep in line with our actual implementation 
m Section El we shall in fact make stronger assumptions, as follows. 

Recah that Q e Ol[X,T] with 2g+l := deg^(Q), that f(r) is the Sylvester 
resultant of Q and ^ w.r.t. X, and r the monic factor of r which is the 
denominator of the connection matrix B{T) = &(r)/r(F) when in lowest terms. 
Define h :— degp(Q). 

We assume that '"(F) modp is squarefree, and that the Laurent expansion 
of B{r) has only negative terms. We say then that S(F) has only simple poles 
modulo p. This ensures that the algorithm in the proof of Theorem 14 . 21 works . 
Let us assume that the local monodromy eigenvalues around each singular point 
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are prepared, so condition (Prep. Rat.) is met and we may apply the precision 
loss bounds in Theorem l4.8l 

To obtain a nice basis for eI'I-, let us further assume that the local mon- 
odromy around the finite poles is nilpotent, and that zero does not occur as a 
local monodromy eigenvalue around the pole at infinity. In this case we may 
take as our basis for E^'l^g the elements {bik} where < i < d—2 {d := deg(r)) 
and 1 < fc < 2g, and the element bik G H^^g{X / S) is the column vector with 
zeros in positions j 7^ k, and in position j — k the 1-form T^dT/r(T). Note that 
the dimension of this space is 2g{d— 1). 

For the complexity analysis, we shall need bounds on the height of the local 
monodromy eigenvalues. Let us assume that a common denominator for the 
local monodromy eigenvalues around each singular point is 2{2g + 1), and when 
written w.r.t. this denominator the numerator does not exceed h(2g — 1) in 
absolute value. Under the assumption that B{T) has only simple poles modulo 
p, we believe that one may prove that the bound on the denominator should 
always holds by a topological argument. Likewise, the author expects that the 
bound on the numerator should also hold, although offers no proof of this. 

Note 7.5 We point out that "generically" in any nice family of polynomials 
both r(r) and f{T) are squarefree and have equal degree. In this case, one ob- 
serves experimentally, and expects to be able to prove, that all residue matrices 
around finite poles are nilpotent. However, the assumption that the degree in 
r of the connection matrix B{r) is less than zero does not hold generically. 
For any family of polynomials Q, e.g. with fixed Newton polytope, one can 
calculate restrictions on the coefficients which must be met. The author has 
no idea of the geometric significance of this assumption. When the assump- 
tion does hold, the local exponents at infinity are observed to exhaust the set 

{ 2(^11 < < 25 odd}. 

7.5 Analysis 

We shall use soft-Oh notation, to hide logarithmic factors in the time and space 
complexity HIl Def. 25.8]. 

Let No be a positive integer which depends upon the equation = Q{X, T) 
in some manner — we shall specify precisely how later. Assume that 

5' (1 + ^ + logpigh)^^ = 0{No), (30) 

i.e., the integer Nq grows at least as fast as the expression on the lefthand-side 
as g,h and p vary. 

7.5.1 Numerical approximations 

Assume that one wishes to compute a p^'' -approximation to the pth power 
Frobenius matrix F : H^igiX) H^^g{X). Then Theorem |01 equation ^ 
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in Note 14.111 inequality (|28|l . Propositions 17.21 and 17.31 and the local mon- 
odromy assumptions in Section 17.41 show that it suffices to take the initial p- 
adic accuracy Nj such that Nj — {a" \ogp{Ni) + P") > Nq for some effective 
constants a" , (3" > 0. For implementations one needs to compute the loss of 
accuracy precisely; however for our complexity estimates it is enough to ob- 
serve a" = 0{g^\ogpig)) and /3" = 0{g^{l + {h/p) + \ogp{ghf)). Here are 
more details. First, Propositions 17.21 and 17.31 combined with the observation 
{6 + 6' — l)degp(f) — 0{gh) shows the following: the F-adic accuracy needed 
in solution of the differential system in Step 1 is 0{pghNj). From inequality 
(|^ . the loss of accuracy in this step is 0{g\ogp{pghN i)) — 0{g{logp{Ni) + 
\ogp{gh))) . Second, the maximum pole order encountered in Step 2 is 0{pghN) 
where N < Nj is the intermediate accuracy, so the loss of accuracy in Step 2 
is 0{alogp{pghN) + /3) = 0{alogp{Nj) + alogp{gh) + f3) where a,f3 are the 
constants in Theorem 14.81 From equation (|21|1 we have a = 0{g^ logp{g)) and 
f3 — 0{g^{l + {h/p) + logp{gh)). Our claim on the loss of accuracy now follows. 
Moreover, from equation H3U|I we see that the initial p-adic accuracy Nj satisfies 
Ni = d{No). 

7.5.2 Time and space complexity 

We now give a precise complexity analysis of the time and space required to 
compute a numerical approximation to the pth power Frobenius matrix F : 
H^^g{X) H^^g{X) using the algorithm in SectionlHI 

Theorem 7.6. Let the affine surface be defined as in Section \7.1\ and assume 
the local monodromy conditions specified in Section \7.4\ hold. We recall that X 
is an open subset of the smooth surface defined by the equation Z'^ = Q{X^ F) 
over the field ¥q of characteristic p, and 2g + 1 := degx{Q), h :— degp((5). 
Let the positive integer Nq satisfy the growth condition Then one may 

compute a p^° -approximation to the pth power Frobenius matrix F : H^^g{X) —> 
H'^^g{X) via the algorithm in Section\^in 0{NQg^h^p\og{q)) hit operations, 
using 0{NQg^hplog{q)) bits of space. 

Proof. Since Nq satisfies growth condition (|30|l . from Scction lV.S.ll we see that 
the initial p-adic accuracy Nj satisfies Nj — 0{No). For the purposes of the 
complexity analysis, we shall forget about the intermediate accuracy N , with 
No < N < Nj mentioned in Section 16.51 and just assume we work with p^'- 
approximations throughout the algorithm. 

Step 1: Using the estimates from Propositions 17.21 and 17.31 we see that the 
F-adic accuracy required in Step 1 is 0{Nipfj,), where 

H := max{degr(g),degr(r),degr(Adj(M))} = Oihg). 

We consider the time/space required to compute an approximation to C{T) in 
Section FS.l.ll The coefficients of F are g x g matrices, whose entries are p^°- 
approximations of elements of the p-adic field K. Moreover, the growth bounds 
given in the analysis following Theorem 15 . II shows that each coefficient requires 
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0{\og{q)Ni) bits of space. This gives a space requirement of 0{Nfg^hp\og{q)) 
bits. For the time, we observe that recurrence 126|l has length bound by 
max{deg(&) + l,deg(r)} ~ 0{gh), and involves multiplication of g x g matrices. 
Thus the time to compute an approximation to C{T) is 0{Njg'^~^'^h^p\og{q)) 
bit operations. One may further compute the approximation to the local Frobe- 
nius matrix F{r) in this time/space, using H23|l . Using the fast radix conversion 
algorithm in Alg 9.14], these time and space estimates are enough for the 
the analytic continuation and radix conversion steps required to make the input 
suitable for Step 2. 

Step 2: The matrix F has size {d - 2)2g = 0{g^h), where d = deg(r(r)). 
Thus 0{g^h) applications of the reduction algorithm from Section ^ are re- 
quired. It is time saving in terms of the parameter g to precompute the inverses 
^\—ir'Im + b)~^, {i'lm + for i and £' in the necessary ranges, as these 

do not depend on the element being reduced. The number of the former in- 
verses is 0{Nip) and each inverse takes 0{Ni \og{q) x ghx g^) bit operations to 
compute; the factor gh arising since the inverse is computed modulo the poly- 
nomial r(r) which has degree 0{gh). There are 0{Nipgh) of the latter inverses 
to compute, but each only requires 0{Nj\og{q)g^) bit operations. Thus pre- 
computation of the matrix inverses takes 0{Nfg'^hp\og{q)) bit operations and 
one needs 0{Nfg^hp\og{q)) bits to store them. The reduction of finite poles 
requires 0{Njp) steps, and each step taking 0{Njg^h\og{q)) bit operations; 
reduction of the pole at infinity requires 0{Nipgh) steps, but each step only 
taking 0{Nig'^ \og{q)) bit operations. The time for the reduction of forms is 
thus 0{Nfg^h'^p\og{q)) bit operations, and this step requires Oh{Nj g^ p\og{q)) 
bits of space. (The time without precomputation of matrix inverses would be 
0{Nfg^h^p\og{q)) bit operations.) This completes the proof. □ 



8 The zeta function of a compact surface 

This section is a direct continuation of Section [3 In particular, throughout 
this section we retain the definitions and assumptions in the preamble to that 
section, as well as those given in Sections 17.11 and 17.41 



8.1 The zeta function of the open surface 

In this section we consider the zeta function Z{X, T) of the smooth afRne surface 
X over ¥q. The trace formula in rigid cohomology for smooth afhnc varieties 
shows 

^^""'^^^ P,{X^T)P,{X,T) 

where P^{X,T) := det(l - Tg^i^^ e 1 + TQp[T]. Certainly 
H'^^g{X) is a one-dimensional Q^-vector space, and Po{X,T) = (1 — q^T). 

Proposition 8.1. Let the polynomial Pi{S,T) be the numerator of the zeta 
function of the open subset S of the projective line; so Pi{S,T) is a product of 
cyclotomic polynomials. Then Pi{X, T) ~ Pi{S, qT) G 1 + rZ[r]. 
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Proof. It is enough to consider the terms i?2.'rig ^^'^ ^^rig spectral se- 

quence for X ^ S, c.f. SectionEHand HTI Eqn (17)]. We have 

:=ker(Vt)-ker(V). 

The isomorphism follows from [31 Cor. 2.6]. We claim the latter space is zero- 
dimensional: Let V € ker(V). Recalling from Section [7^ that the local mon- 
odromy eigenvalues around finite poles are all zero, expanding v around the 
finite poles one deduces that v G K^^ . Since the local monodromy eigenval- 
ues around the pole at infinity are non-zero, expanding v around this pole one 
deduces v = 0. We have 

■■= coker (^A ; hI!,^{X/S) ^ H%{X/S)dr^ . 

But H^^g{X / S) = A\ the weak completion of the coordinate ring of S. So 
det(l - Tg2F-i°g.(«)|i?i;°^^) = det(l - (Tg)gF- = P^{S,qT). 

□ 

Proposition 8.2. The polynomial P2{X,T) has integer coefficients. 

Proof. Integrality follows from Proposition 18. II since the zeta function itself is a 
power series with integer coefhcients. □ 

Kedlaya's p-adic analogue of Deligne's main theorem tells us that the com- 
plex absolute values of reciprocal zeros of P2{X, T) belong to the set {1, q^'"^, q} 
|2(J| : we will deduce this is an elementary manner in Section [8.2l 

8.2 The zeta function of a compactification 

In this section we show that the full zeta function Z{X, T) of a compactification 
X oi X may be easily recovered given the first 0{gh) coefficients of P2{X,T) 
to precision modulo where N = 0{gh\og{q)). 

To simplify the analysis, and keep in line with our actual implementation, 
we shall make some further restrictions on the polynomial Q{X, F). Recall that 
Q is monic in X of degree 2g + I, and has degree h inV. Let us further assume 
that it is monic in F with h odd, has constant term 1, and that 2g + 1, h and 
the prime p are mutually coprime. Moreover, assume that all other terms in 

— Q{X, F) have exponents lying within or on the boundary of the polytope 
A with vertices the origin and the points {2g + 1, 0, 0), (0, ft., 0) and (0,0,2). 
Then the Newton polytope of — Q{X, F) (taken modulo any prime number 
p) is the simplex A. We assume that Z^ — Q{X, F) is non-degenerate w.r.t. 
the polytope A c.f. |3 Sec. 3.6]: Specifically, the polynomials Q{X,0),Q{0,T) 
are squarefree, and Q, ^ have no common solutions. Let X be the toric 
compactification of the affine variety X in the toric projective space Pa c.f. 
13 Sec. 3.2]. This is a smooth compact variety. Since the outer face of A 
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is a triangle with no interior points, it follows that X = Xaff U where 
Xaff := Spec(F,[X, T, Z]/iZ^ - T))). 

One does not need to be familiar with the exact details of the construction: 
the point is simply that we have compactified the zero set in afline space of the 
equation = Q{X,r) by adding a single projective line. 

Definition 8.3. Let P{T) e Z[r], q a prime power, and £ a non-negative inte- 
ger. We call P{T) pure of weight i with respect to q if its reciprocal zeros have 
complex absolute value q^l"^ . We shall just say P is a weight i Weil polynomial, 
when q is understood. 

Proposition 8.4. Let X be the smooth toric compactification of the set of affine 
solutions of the equation Z^ — Q{X,Y), as defined immediately above. Then 
the zeta function of X has the form 

^^^'^^ " {l-T)P2{X,T){l-q^T)' 
where P2{X,T) G Z[T] is a Weil polynomial w.r.t. q of weight 2, and 

deg(P2(^,T)) = r(2A) - 4/*(A) - 3 - ^(r(A') - 1). (32) 

A' 

Here the function I* counts lattice points in the interior of a polytope or polygon, 
and the .sum is over the two-dimensional faces A' of A. 

Proof. The claim on the weight follows from Deligne's theorem [Hj. The other 
claims follow from the formula for Hodge-Deligne numbers of complex toric sur- 
faces in 13 Sec. 5.11(c)], the comparison theorem between singular cohomology 
and Z-adic etale cohomology for the modular reduction of smooth complete va- 
rieties over number fields, and the trace formula in £-adic etale cohomology. □ 

We note, but do not use, that the formula for deg{P2{X ,T)) is valid for 
arbitrary Newton polytopes A, assuming that Pa is smooth and Z^ — Q{X, T) 
is non-degenerate w.r.t. A. 

To simplify the statement and proof of the next theorem, and again to keep 
in line with our implementation, we make some further restrictions: Assume 
that the Sylvester resultant f{T) is squarefree modulo p of degree d = degp(r), 
and that for each 7 G with f{'~f) = mod p, the "missing fibre at F = 7" in 
the pencil X ^ S has a unique double point. 

Proposition 8.5. Definitions and assumptions as in Sections \ 7.1\ \7.4[ 
and the present section. Let X = X U C where C is a union of curves. We 
recall that X is a compactification of the smooth surface defined by the equation 
Z^ = Q{X,T) over ¥g, with 2g -f 1 := degx(Q) and h := degr(Q). Th 
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Z[C,T) may be computed deterministically in O{g^hp\og{q)^) bit operations. 
Moreover, given Z{C,T) the zeta function Z{X,T) may be recovered from the 
first d — 2g coefficients in P2{X, T) each to p-adic precision modulo p^ where 

max (log (2(7' ('^ . 



N := 
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Proof. Define f{T) := r{T) mod p. Let r(r) = 11^=1 ^i(r) be the irreducible 
factorisation and define di :— deg(fi). For i = 1, . . . , s denote 7i := F G Ki := 
Fg[F]/fi(F). Our assumption tliat each singular fibre has a unique double point 
implies that Q{X,jt) = [X ~ ai)'^Hi{X) where Hi{ai) ^ 0. Define Si 



if Hi{ai) is a square in Ki, and 5i 
follows that 

Z{C,T) 



-1 otherwise. Since X = Xaff U '. 



-1 
it 



n 



(33) 



where Pi{T) is the numerator of the zcta function of the genus .9 — 1 curve 
Z'^ = Hi{X). Note that Z{X,T) = Z{X,T)Z{C,T). From (EU, Propositions 
18.11 and 18.41 and (|33|) . and by noting the weights of the different factors, we 
deduce the following: 



P2iX,T) = 
P2iX,T) = 



W2{P2iX, T)) n:^i RdT^^){l + S,T^'') 

{l-qT)w2{P2{X,T)). 



Here W2{P2{X^T)) is the "interesting" weight two factor in P2{X,T). It has 
degree 2g{d — 1) — (d + d{2g — 2)) — d — 2g. The theorem now follows, using 
Kedlaya's algorithm JS] to compute Z{C, T) and noting d = 0{gh). □ 

Note that W2{P2{X,T)) satisfies the same functional equation as P2{X,T). 
The sign in this functional equation is (—1)* where s is the multiplicity of {1+qt) 
as a factor of P2{X, T). This is unknown. However, by computing only the first 
[(d — 2g)/2\ + 1 coefficients in P2{X, T) to p-adic precision modulo where 



N 



logj, 2q 



d-2g 



, e := 



d-2g 



(34) 



one can find two possible candidates for P2{X,T). One hopes that only one is 
a weight 2 Weil polynomial! 



8.3 Computation of the full zeta function 

In this section we retain the definitions and assumptions in Sections 17.11 17.41 
and Section Theorem 17.51 and Proposition 18 . 51 together yield an algorithm 
for computing the full zeta function of the compact surface X, provided we can 
estimate the loss of precision between the computation of the absolute Frobe- 
nius matrix F : H^^g{X) H^^g{X) and the calculation of coefficients in the 

polynomial P2{X,T) = det(l - Tq^F-^°iip'^'}^\H^^giX)). Note that in practice 
one actually computes coefficients in the polynomial det(T — i^'°Sp('?) \H^ig{X)). 

Theorem 8.6. Fix a positive constant C and positive integer g. Assume that 
degx{QiX,T)) = 2g + l and that h := d_egr(Q(X, F)) satisfies h/p < C. Then 
one may compute the zeta Junction Z{X, T) of the compactification X of the 
affine surface defined by Z^ = Q{X,T) in (D{h'^plog{q)^) bit operations using 
0{h^p\og{q)^) bits of space. 
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Note that the hidden constants in the Soft-Oh notation depend upon both 
the genus g and constant C. 

Proof. Since g is fixed and h/p is bounded, the numbers a and /3 in Theo- 
rem 14.81 a' in Theorem 15.11 and consequently a" and /3" in Section 17.5.11 are 
bounded absolutely, independent oi X. It follows easily from Theorem 14.81 and 
Propositions 17.21 and 17.31 that the Frobenius matrix F has valuation bounded 
below by some absolute constant — c, with c > 0. We require the final p-adic 
precision to be modulo with N as in the statement of Proposition l8.5l Notice 
= 0{\og{q)gh). A naive analysis of the loss of accuracy during the computa- 
tion of the characteristic polynomial from the absolute Frobenius matrix shows 
that it certainly suffices to take Nq — N + c\ogp{q){d — 2g) + ordp((d — 217)!) 
in Theorem 17.61 recall d = 0{gh). Note that condition (|30|l is trivially sat- 
isfied in this case since the LHS is bounded absolutely. The complexity es- 
timate follows by putting this value for Nq in Theorem 17.61 and noting that 
the resulting time/space estimate also suffices for computing the characteristic 
polynomial. □ 

We note that assuming Conjecture 17.41 is true, and that the pole order is 
bounded in some manner depending only on g, them Theorem 18 . 61 holds without 
the restrictions on the relative growth of h and p. The point is that in this case 
one can take a" and (3" to depend only on by using the original Christol- 
Dwork theorem, see Note 14. ill 

The author has been unable to prove that the valuation of the Frobenius 
matrix F is bounded below by some absolute constant, i.e., bounded indepen- 
dent of g and h. If one could show this then putting Nq = O{gh\og{q)) in 
Theorem 17.61 and assuming h/p remains bounded as h and p vary, we get the 
estimate 0{g'^ h'*^p\og{q)^) bit operations/(!)(5^/i'^plog(q)'^) bits of space, for the 
computation of the zeta function Z{X,T). Note that since the middle betti 
number d ~ 2g + 1 in this case is approximately gh, this compares favourably 
with the "deformation algorithm"; see the end of Section 13 

9 Surfaces: implementation and experiments 

In this section we report on a Magma (v. 2. 11-2) implementation of our algorithm 
for the surfaces described in Sections [7| and |S1 in the case of a prime field. The 
experiments detailed were peformed using a dual processor Intel Pentium 4 
(3GHz with 1024Mb cache and 2Gb RAM per processor). Time and space 
requirements stated are as returned by the in-built Magma function. 

9.1 Examples 

All of the examples satisfied the hypothesis in the statement of Proposition l8.5l 
and the local monodromy eigenvalues were as observed in Note 17.51 
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Example 9.1 Let 



Q{x, r) := X'''+(4r''+5r''')x+r"+6ri2+5ri°+8r^+8r^+5r^+r'*+5r''^+r2+i, 

and p := 17. Then the Sylvester resultant r{T) := Res(X, Q. is squarefree 
modulo p and equals, up to a constant, the denominator r{T) of the connection 
matrix 6(r)/r(r). Both polynomials have degree d := 2 x 13 = 26. The genus 
of the generic fibre is g := 1. The space H^^g{X) associated to the open surface 
^2 =_Q(X,r),r(r) ^ Omodp has dimension {d - l)2g = 50. The space 
H^^g{X) associated with the smooth toric compactification has dimension d — 
2g+l = 25. Computing a matrix for the Frobenius map F : H^^g{X) H^^g{X) 
to precision modulo p^^, we recovered two possible choices for the polynomial 
det{T — F\H^^g{X)). Only one was the reciprocal of a weight 2 Weil polynomial 
w.r.t. 17. Specifically, P2{X,T) = (1 - 17r)2(l + 17T)R(T) where R{T) is the 
irreducible polynomial 

1 + 2^r^ + 2^3^17^r^ + 2*17h9^r^ + 2^3^5h7V + 3h7^23^r^ + 17*23^7" 

-2^5^ 17^r'^ + 3''l7®r* + 2^5^7^ 17'^r^ + 2^ 17** 191^r"' + 2*13^ 179^11 + 2^ 17^" 191^r^^ 
+2^5^7h7"r^'V3''l7^^r"-2^5h7"r^^+17^'*23^r^V3h7^^23^r^V2^3^5h7^'^r" 

The Hodge numbers defined a polygon called the Hodge polygon which lies 
below the Newton polygon of P2{X,T). In this case, the Hodge numbers are 
2,21,2 which explains the high divisilbility of the coefficients by powers of p, 
c.f. P Remark 1.6.4]. 

The computation is provably correct under no additional hypothesis. It 
took just under 23 hours and 13 minutes, and required just under 1.312 Gbytes 
of memory. We note that over half the time required was taken computing 
r(r)-expansions of the elements in the our relative Frobenius matrix. This was 
necessary to ensure the input to the second stage of the algorithm was in the 
appropriate form. 

Under Conjecture 17.41 this example required just under 2 hours 36 minutes, 
and 216Mbytes of memory — the pole order appears to be 39. 

Example 9.2 Let 

Q{x, r) := x^ + (r" + + i)x + + 2T^^ + + + 2r + 1, 

and p := 5. So rf = 2 x 31 = 62, g := 1, dim{H^,g{X)) = 2g{d - 1) = 121, 
dim{H^^g{X)) ~ d — 2g + 1 = 61. The characteristic polynomial dct(T — 
F\H^^g{X)) was computed modulo p^^. We found P2{X,T) ^ {1 - 5T)'^{1 + 
5T)R{T) where R{T) is the irreducible integer polynomial 

1 - 5^T^ + 2^3^5^T^ - 2^3^5^r^ - 2^5^r* + 5^89^r^ - 2^5^409^7^ + 2^5^7U7^r'^ 

-3^5^617^r* - 5*37^r^ + 5^7^727^r^° - 3^5^11^97^r" + 5*11^53^r^^ + 5'-*7^73^r^^ 
-5''7^43^r^'* + 5"ll^61^r^^ - 2^5"829^T^'' - 2^5^^677^T^'^ + 2^3^5"79^T^* 
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+2^3^5^°29^r^^ + 5^^367^r^'' - 5^^7^23^r^^ + 2^3^5^^53^r^** - 2^3^5^''7n3^r^® 
+2^3^5^^53^r^°-5^^7^23^r^^+5^^367^r^^+2^3^5^*29^r^^-3^5^''7^r^''+5^^137^r^'^ 
+2^3^ 5^^ 11^ T^^ - 5^''37^r^^ + 2^5^^1777^r^** - 5^*'53^89^T-''-' + 2^3^5^®79^r"*° 
-2^5^«677^r*^ - 2^5''829^r*2 + S'^ll^Gl^T^^ - 5^^7'43^r^^ + 5*^7^73^r^^ 
+5"'ll^53^r"« - 3^5*'ll^97^r*" + 5"='7^727^r"» - 5^^37'r*^ - ^'b^^n^T"" 
+2^5'*^7M7^r^^ - 2^5''*409^r^^ + 5™89^r^^ - 2^5^^r^* - 2^3^5^^r^^ + 2^3^5^^r*'^ 

The Hodge numbers in this case are 5,51,5. 

The computation is provably correct only under Conjecture 17.41 — the pole 
order appears to be 31. It took 13 hours and 577 seconds, and required just 
under 834 Mbytes of memory. 



Example 9.3 Let 

Q{x, r) := + ix^ + (4r2 + 4r + 8)x + + 5r^ + 1 

and p 11. So d ^ 4 x 7 = 28, (/ 2, dim{H^^g{X)) = 2g{d - 1) = 108, 
dim{H^^g{X)) ~ d — 2g + 1 = 25. The characteristic polynomial det(T — 
F\H^.^g{X)) was computed modulo p^^. We found P2{X,T) = (1 - llTfRiT) 
where R{T) is the irreducible integer polynomial 

1 + i9^r^ + i7^r^ + 2-'"'ii^r^ + 7^ii^r'* + 3Hi^i03^r^ + 2^3hi''4i^r^ + 3^ii'^23^r^ 

+ 11® ISl^T* + 2-'^5^11^r'' + 2^11** 47^T^° + 3Hi^59^T" + 2^11^'' 47^r^^ + 2^5^1l"r^^ 

+ii^^i5i^r" + 3^ii^^23^r^-'^ + 2^3ni"4i^r^'' + 3^ii^^i03^r^'^ + r^uisyis 

j q5 -j^ 18 j -j^ -j^ 18 )2-i20 j -j^ 2 g 1 ^2"' ^ | -|^-j^22^2~'22 

The Hodge numbers are 2,21,2 in this case. 

The computation is provably correct only under Conjecture 17.41 — the pole 
order appears to be 21. It took just under 14 hours 36 minutes and required 
4.41Gbytes of memory. 

We note that use of Method 2 fSection l5.1.2|l rather than Method 1 (Section 
I5.1.1|l significantly reduces the space requirement; however, we do not have 
provable precision loss bounds for Method 2. If one is satisfied with plausible 
rather than provable output, larger examples may be computed. 



9.2 Calculation of precisions required 

We now address the delicate problem of minimizing the amount of precision 
one needs to carry through the algorithm to obtain an answer which is provable 
correct (possibly assuming Conjecture 17. 4|) . 

Fix a positive integer and suppose that we wish to compute a p^^- 
approximation to a matrix for F acting on H^^g{X). We will discuss the choice 
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of A^3 later in this section. Let i?2g,p be as in H2U|I : in particular, for p > 2g we 
have i?2ff,p = 2(7 — 1. Recall that 2g + 1 := degxiQ)- 

Define Xfm and N2jin to be the smallest integer solutions to the inequalities: 

lXfin/p\ - l\ogpi2Xfin + l)\ > N2Jin 

iV2,/^„- (2B2g,p + 2g)[logp(a;/„0J > N3. 

More precisely, let x fin be the smallest integer solution to 

[xf^n/p\ - Llogp(2a:/„, + 1)J - (2S2g,p + 25)[logp(a;/,„)J > N3, 

and define A^2,/m in the obvious way. Applying Proposition 17.21 and Theorem 
14.81 and recalling that we have nilpotent monodromy around the roots of r{T), 
one sees the following: It is enough to compute the coefficients /^(r) with fc < 
in the r(r)-adic expansion of the entries of i^(r) for |fc| < Xfm, and to compute 
these with p-adic precision "modulo p^^.f^^" . The point is that for any basis 
form bik(r) fSection l7.4|) . the coefficients in "(i, A;)th column" of the matrix for F 
are given by applying the reduction algorithm from the proof of Theorem l4.2l to 
the image form F{T)bik{TP)', but the reduced form to p-adic precision "modulo 
is not affected by coefficients /fc(r) for k negative with |fc| > Xfin- Define 
Nrjin ■= deg(r)x/j„. 

One can argue in a similar manner to determine which coefficients /fc(r) 
for fc > in the r(r)-adic expansions of entries in F{T) must be computed, 
and to what precision. We return to this shortly, but let us say that we have 
determined suitable integers Xinf and N2,inf, and defined A^r,inf := deg(r)a;m/- 

Define N2 := max{iV2,/m, A'2,in/}. We need to compute the coefficients 
/fe(r) in the r-adic expansion of entries in the global Frobenius matrix F(T) 
for —Xfin < k < Xinf with p-adic precision "modulo p^^" . Define Nr ■= 
Nrjin + Nr^inf- Since there is no loss of accuracy during the analytic continu- 
ation stage f Section 15. 2|) . it is enough to compute a -approximation to the 
local Frobenius matrix F{T) modulo F^^^ Using the method in Section [5.1. 11 
equation (|28|l tells us we must perform the local calculation itself to p-adic 
precision "modulo p^^ " , where 

iVi := N2 + i3B2g,p + l)Llogp(iVr)J - B2g,p + min{ordp(F(0)), 0}. 

Note that ordp(F(0)) > when p > 2g. Our algorithm begins by computing 
a p-'^i -approximation to the matrix F(0); see ^H] for an analysis of the loss of 
accuracy during this initial computation. 

We return to the question of determining N2^inf, a;inf and N-p^inf- One can 
do this via an analogous system of inequalities to those above, using Proposition 
17.31 and Theorem 14.81 combined with the general estimates for a and /3 derived 
from Note 14.111 The problem is that since the local monodromy around the 
point at infinity is not nilpotent, the constants a and (3 are rather large. The 
get around this, the author wrote a short computer program which calculated 
more careful bounds on the growth of the coefficients in the uniform part of the 
local solution matrix around the point at infinity. In the notation of Lemma l4.9l 
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2x,r,f-{p-l){S-l) 

2p{5+S'-l) 



the author computed a lower convex function ai{i) such that ovdpiYi) > —ai{i) 
for all J > 1. The function ai(i) depended explicitly on the local monodromy 
eigenvalues at infinity and p; the time required to compute ai(i) grew as logp(z) 
with i > 1. Here are brief details: For eigenvalues in the interval [0, 1) use 
[TUl Lines 8-9,19, Page 196]; for general prepared eigenvalues, use the proof of 
Lemma f4. 91 but compute a tighter lower bound on "ordj,('hI)" via the proof of 
Lemma 14.101 and use the inequality "ordp(Fi) > ordp(yi+2A) + ordp(7Y)". The 
function ai{i) was fed as input to the analysis in the proof of Theorem 14.81 to 
yield a better function a{£), say, which could be used on the righthand-side on 
the statement of the theorem. With this more refined function, one takes Xinf 
and N2,inf to be the smallest integer solutions to the inequalities: 

{Vtnf + 1) - Ll0gp(p(2yj„/ + 1))J > N2anf-, {Vtnf 

deg(r)) > N3. 

See (|29(l for the definitions of the numbers S and S' . 

When the author assumed Conjecture 17.41 was true, he did not perform the 
calculation in the preceding paragraph, but instead defined iV2,,;„/ :— N2,fin 
and iVr,j«/ 100. 

We require that the characteristic polynomial dct(r — F\H'^.^g{X)) be com- 
puted modulo p^, with N as in equation (|34|) . If ordp(F) > 0, then one can 
take := TV; this was the case in Example O If ordp(i^) < 0, there may 
be some loss of accuracy during the computation of the characteristic polyno- 
mial. The author had an ad hoc solution to this problem: Specifically, it was 
observed in practice that even when ordp(_F) < 0, some small power of F had 
non-negative or even positive valuation. By examining the valuation of powers 
of F, and using the formula P/^{X,T)/ P2{X,T) = - Tr(F'=)r'="i, one 

can deduce explicit bounds on the loss of precision. This enabled the author 
to establish usable and provable precision loss bounds during the calculation of 
the characteristic polynomial; however, when the initial computation revealed 
ordp(i^) < 0, one did need to rerun the computation with an increased value for 
A'3 to get a provably correct answer. 

The parameters [N, N3, N2jin, N2,inf, Ni; N-pjin, Nr.inf] in the examples 
were set as follows: in Example EH [18, 18, 26, 56, 67; 12376, 16692] uncondi- 
tionally and [18, 18, 26, 26, 37; 12376, 100] under ConjectureEl in Example lO 
[55, 60, 72, 72, 95; 23560, 100]; in Example ESI [19, 25, 45, 45, 72; 14476, 100]. 
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